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20  Key  Points 

21  •  An  RMM-like  index  was  created  for  the  skeleton  model  that  mimics  observations 

22  •  Stochasticity  helps  improve  MJO  initiation  and  termination  event  statistics 

23  •  The  skeleton  model  produces  more  realistic  MJOs  when  forced  with  observed  SSTs 

24 
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25  Abstract 

26  The  Madden-Julian  oscillation  (MJO)  skeleton  model  is  a  low-order  dynamic  model  that 

27  is  capable  of  simulating  many  of  the  observed  features  of  the  MJO.  This  study  develops  a  model- 

28  based  “MJO”  index  that  is  similar  to  the  well-known  real-time  multivariate  MJO  (RMM)  index 

29  to  better  facilitate  comparison  between  the  skeleton  model  and  observational  data.  Multivariate 

30  and  univariate  empirical  orthogonal  function  (EOF)  analyses  were  performed  on  the  convective 

3 1  heating  and  zonal  wind  data  taken  from  the  skeleton  model  for  simulations  forced  with  an 

32  idealized  warm  pool  and  observed  sea  surface  temperatures  (SSTs).  The  leading  EOF  modes 

33  indicated  a  wavenumber- 1  convectively-coupled  circulation  anomaly  with  zonal  asymmetries 

34  that  closely  resembled  the  observed  RMM  EOFs,  especially  when  the  model  was  forced  with 

35  observed  SSTs.  The  RMM -like  index  was  used  to  compute  an  MJO  climatology  and  document 

36  the  occurrence  of  primary,  continuing,  and  terminating  MJO  events  in  the  skeleton  model.  The 

37  overall  amount  of  MJO  activity  and  event  lengths  compared  reasonably  well  to  observations  for 

38  such  a  simple  model.  Attempts  at  reconciling  the  observed  geographic  distribution  of  individual 

39  MJO  initiation  and  termination  events  were  not  successful  for  the  stochastic  simulations,  though 

40  stochasticity  is  necessary  in  order  to  produce  composite  MJOs  that  initiate  and  decay  with  time 

41  scales  similar  to  observations.  Finally,  analysis  indicates  that  the  existence  of  slow  moving, 

42  eastward  travelling  waves  with  higher  wavenumbers  (k  ~  12)  embedded  within  the  large-scale 

43  flow  often  precedes  MJO  termination  in  the  skeleton  model. 
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44  1.  Introduction 

45  Intraseasonal  variability  of  the  tropical  atmosphere  is  dominated  by  an  eastward 

46  propagating  equatorial  mode  known  as  the  Madden- Julian  oscillation  (MJO),  appropriately 

47  named  after  its  initial  discoverers  [ Madden  and  Julian ,  1971;  1972].  The  MJO  represents  a 

48  convectively-coupled  mode  of  cloud  and  circulation  anomalies  on  the  planetary  scale  (i.e., 

49  wavenumbers  1-3)  with  an  oscillation  frequency  of  roughly  30-60  days  [ Madden  and  Julian, 

50  1994;  Zhang,  2005].  Deep  convection  associated  with  an  MJO  event  generally  appears  as  a  broad 

5 1  envelope  of  cloud  and  precipitation  anomalies  that  most  commonly  begins  over  the  western  or 

52  central  Indian  Ocean  [e.g.,  Matthews,  2008;  Straub,  2013]  and  propagates  eastward  at  ~5  m  s"1 

53  until  eventually  terminating  over  the  maritime  continent  or  the  central  Pacific  [e.g.,  Stachnik  et 

54  al.,  2015].  The  MJO  directly  impacts  local  weather  and  tropical  climate  and  has  also  been  shown 

55  to  have  far  reaching  dynamical  effects  including  modifying  the  global  monsoons  [Lau  and  Chan, 

56  1986;  Goswami  and  Ajaya  Mohan,  2001;  Matthews,  2004;  Wheeler  and  Hendon,  2004],  genesis 

57  and  intensification  of  tropical  cyclones  [Maloney  and  Hartmann,  2000;  Jiang  et  al.,  2012],  mid- 

58  latitude  weather  anomalies  [Higgins  et  al.,  2000;  Bond  and  Vecchi,  2003;  Guan  et  al,  2012],  and 

59  may  impact  the  onset  of  El  Nino-Southern  Oscillation  (ENSO)  events  [e.g.,  Lau  and  Waliser, 

60  2005], 

61  Despite  its  global  importance,  successfully  simulating  the  MJO  remains  a  challenge  in 

62  most  global  climate  and  weather  forecast  models  [e.g.,  Slingo  et  al,  1996,  Waliser  et  al.,  1999, 

63  Zhang  et  al.,  2006;  Hung  et  al.,  2013;  Jiang  et  al.,  2015].  Global  climate  models  (GCMs) 

64  generally  struggle  with  initiating  MJO  events  and  often  simulate  weaker  eastward  propagation 

65  than  observed  in  nature.  Models  with  more  realistic  MJOs  usually  tune  specific  model 

66  parameters  including  estimates  of  air-sea  coupling,  the  sensitivity  of  parameterized  convection  to 
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67  atmospheric  moisture,  or  increase  deep  convection  entrainment  rates  in  order  to  delay  convection 

68  (among  others),  which  often  comes  at  the  sacrifice  of  a  larger  bias  in  the  mean  state  [e.g., 

69  Bechtold  et  al.,  2008;  Hannah  and  Maloney ,  2011;  Benedict  el  al.,  2013;  Klingaman  and 

70  Woolnough,  2014a;  Klingaman  and  Woolnough,  2014b].  Recent  advances  using 

7 1  superparameterizations,  in  which  a  cloud  resolving  model  is  embedded  within  each  GCM 

72  gridpoint,  have  demonstrated  some  success  at  improving  simulations  of  the  MJO  [e.g., 

73  Grabowski ,  2003;  Benedict  and  Randall ,  2009;  Pritchard  and  Bretherton,  2014].  Nevertheless, 

74  these  enhancements  are  inherently  linked  to  a  significant  increase  in  the  computational  expense. 

75  Simple  theoretical  models,  such  as  the  multicloud  model  of  Khouider  and  Majcla  [2006,  2007], 

76  have  also  been  shown  to  improve  simulated  MJO  variability  when  incorporated  into  aquaplanet 

77  and  coarse-resolution  GCMs  [ Khouider  el  cd. ,2011;  Ajayamohan  et  al. ,  2013;  Deng  et  cd. ,  2015] . 

78  Regardless,  there  is  a  pressing  need  for  the  refined  development  of  computationally  inexpensive, 

79  low-order  models  that  are  fully  capable  of  simulating  the  large-scale,  observed  features  of  the 

80  MJO.  Furthermore,  incorporating  these  fundamental  MJO  dynamics  into  higher-order  models 

81  may  also  help  improve  intraseasonal  and  medium-range  weather  forecasts. 

82  The  MJO  skeleton  model,  developed  by  Majcla  and  Stechmann  [2009,  201 1],  is  one  such 

83  example  of  a  low-order  dynamical  model  that  successfully  simulates  a  number  of  salient  features 

84  of  the  MJO  including  (1)  an  eastward  propagation  at  ~5  m  s"1,  (2)  a  unique  dispersion  relation 

85  with  dco/dk  ~  0,  and  (3)  a  horizontal  quadrupole  vortex  structure.  Recent  work  by  Thual  et  cd. 

86  [2014]  incorporated  a  stochastic  element  into  the  convective  parameterization  of  the  skeleton 

87  model  and  was  able  to  further  simulate  (4)  observed  MJO  event  intermittency  and  (5)  the 

88  organization  of  MJOs  into  circumnavigating  events  and  wave  trains  as  seen  in  nature  [e.g., 

89  Matthews ,  2008;  Yoneyama  et  cd.,  2013].  These  MJO-like  waves  have  also  been  positively 
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90  identified  in  observational  data  by  Stechmann  and  Majda  [2015]  and  Ogrosky  and  Stechmann 

91  [2015]  using  a  new  index  that  projects  reanalysis  data  onto  the  theoretical  linear  solutions  of  the 

92  MJO  obtained  from  the  skeleton  model  (i.e.,  the  amplitudes  of  the  Kelvin  and  equatorial  Rossby 

93  wave  structures,  moisture  anomalies,  and  amplitude  of  convective  heating).  While  the  above 

94  index  is  proficient  at  isolating  the  observed  MJO  signal  without  the  use  of  time  filtering  or 

95  empirical  orthogonal  functions  (EOFs),  its  use  should  be  supplemented  by  more  common  MJO 

96  metrics  in  order  to  better  facilitate  comparisons  with  higher-order  GCMs  and  operational  weather 

97  forecast  models. 

98  This  work  intends  to  evaluate  and  build  a  complementary  index  for  the  MJO  skeleton 

99  model  based  on  the  widely-used  real-time  multivariate  MJO  (RMM)  index  of  Wheeler  and 

100  Hendon  [2004].  The  RMM  index  has  experienced  tremendous  popularity  over  the  last  decade 

101  and  there  are  a  large  number  of  studies  compositing  MJO  data  by  RMM  phase  that  permit  a  more 

102  direct  comparison  between  the  skeleton  model  and  observations  or  GCMs.  Moreover,  defining 

103  an  RMM-like  index  for  the  skeleton  model  provides  the  opportunity  to  evaluate  higher-order 

104  features  of  the  model  solutions  in  a  framework  consistent  with  observations  including  MJO  event 

105  climatology  and  the  precursor  conditions  associated  with  the  initiation  and  termination  of  MJO 

106  events. 

107  This  article  is  organized  as  follows.  Section  2  gives  a  brief  overview  of  the  skeleton 

108  model  in  addition  to  the  methods  used  to  derive  the  RMM-like  index.  We  extend  the  recent  work 

109  by  Ogrosky  and  Stechmann  [2015]  and  run  the  model  using  either  an  idealized  warm  pool  or  a 

1 10  more  realistic  background  state  based  on  observed  sea  surface  temperatures  (SSTs). 

111  Comparisons  of  the  fundamental  MJO  structures  in  the  model  are  evaluated  using  the  RMM-like 

112  framework  and  presented  in  section  3.  We  also  evaluate  the  overall  level  of  MJO  activity  in  the 
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113  model  along  with  several  applications  that  include  a  comparison  to  the  observed  MJO  event 

1 14  climatology  and  identification  of  precursor  conditions  potentially  responsible  for  MJO 

115  termination  following  Stachnik  et  al.  [2015].  Finally,  section  4  closes  with  a  summary  and 

116  discussion  of  key  results  while  identifying  avenues  for  continued  and  potential  future  work. 

117 

118  2.  Data  and  Methods 

119  2.1.  The  Skeleton  Model 

120  The  MJO  skeleton  model  [Majda  and  Stechmann,  2009,  2011]  is  a  low-order,  nonlinear 

121  oscillator  designed  to  qualitatively  reproduce  the  main  features  of  tropical  atmospheric  dynamics 

122  at  intraseasonal  time  and  planetary  scales.  The  skeleton  model  is  unique  in  that  the  large-scale 

123  waves  are  assumed  to  arise  through  neutrally  stable  interactions  involving  convectively-coupled 

124  waves  at  the  synoptic  and  mesoscale  that  are  driven  by  planetary  scale  moisture  anomalies 

125  instead  of  instabilities  that  occur  on  the  planetary  scale.  Strictly  speaking,  the  basic  model 

126  premise  assumes  that  the  instabilities  and  damping  take  place  on  the  synoptic  or  mesoscale,  the 

127  details  of  which  are  unimportant  for  the  resulting  convective  heating  that  drives  the  large-scale 

128  circulation  anomalies  at  the  planetary  scale. 

129  Consequently,  the  skeleton  model  comprises  a  simple  Matsuno-Gill  type  model  of  the 

130  equatorial  primitive  shallow  water  equations,  without  damping,  to  sufficiently  represent  the 

131  planetary  scale,  dry  dynamics  (Equations  1-4).  In  addition,  there  are  two  evolution  equations  that 

132  describe  the  time  tendencies  for  the  amplitude  of  the  convective  heating  wave  envelope  and 

133  lower  tropospheric  moisture  anomalies,  respectively  (Equations  5-6).  Given  that  observations  of 

134  the  MJO  often  demonstrate  positive  low-level  moisture  anomalies  preceding  and  in  quadrature  to 

135  regions  of  deep  convection  [e.g.,  Hendon  and  Salby,  1994;  Myers  and  Waliser,  2003],  the 


7 


136  skeleton  model  parameterizes  the  time  tendency  of  the  convective  envelope  amplitude  (a)  as  a 

137  function  of  lower  tropospheric  moisture  anomalies  ( q ).  This  relationship  is  signified  by  the 

138  nonlinear  amplitude  equation  in  the  following  set  (Equation  7).  The  model  equations  in  non- 

139  dimension  units  are  thus 


140 

du  dp 

dt  gx 

(1) 

141 

dp 

yu  =  -^ 

(2) 

142 

o  =  -^+6 

dz 

(3) 

143 

du  dv  dw 

dx  dy  dz 

(4) 

144 

—  +  w  —  Ha  —  S9 
dt 

(5) 

145 

—  Qw  —  —Ha  +  sq 

dt 

(6) 

146 

da  „ 

«  =  r«a- 

(7) 

147  where  u,  v,  and  w  represent  the  zonal,  meridional,  and  vertical  velocities,  respectively,  p  is 

148  pressure,  and  Oh  potential  temperature.  H  and  Q  are  constants  that  scale  the  heating  rate  and 

149  approximate  the  background  vertical  moisture  gradient,  respectively.  The  quantity  /"is  a  positive 

150  proportionality  constant.  Finally,  sq  and  s°  represent  external  radiative  cooling  and  moistening 

151  terms  that  can  be  used  to  force  the  model.  It  should  be  noted  that  all  variables,  with  the  exception 

152  of  a ,  are  anomalies  from  a  radiative-convective  equilibrium  state. 

153  We  use  the  simplest  version  of  the  skeleton  model  herein,  where  the  above  equations  are 

154  reduced  to  a  single  dimension  by  horizontally  and  vertically  truncating  to  retain  only  the  first 

155  meridional  Hermite  function  and  vertical  baroclinic  mode.  The  effect  of  the  meridional 

156  truncation  is  to  maintain  an  equatorially  symmetric  heating  structure  that  excites  only  Kelvin 
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157  waves  and  the  first  equatorially  symmetric  Rossby  waves.  Details  of  the  truncations,  along  with 

158  the  equations  for  the  model  solutions  and  their  reconstructed  variables,  are  provided  in  Majcla 

159  and  Stechmann  [2009,  2011]. 

160  Recently,  Thual  et  al.  [2014]  modified  the  nonlinear  term  qa  in  the  convective 

161  parameterization  (Equation  7)  to  include  a  stochastic  birth-death  process  to  better  account  for  the 

162  effects  of  synoptic  and  mesoscale  processes  that  introduce  significant  irregularity  to  the  MJO. 

163  Specifically,  the  amplitude  of  the  convective  envelope  was  replaced  by  a  random  variable,  a,  that 

164  evolves  according  to 

165  a  —  A  aq,  (8) 

166  where  A  a  is  a  fixed  step  and  //  is  a  non-negative  integer.  The  probability  (P)  of  any  given  value  // 

167  is  prescribed  by  the  master  equation 

168  =  ^  "  ^P(J1  -V-  +  [b0?  +  1)^0?  +  1)  -  b(??)E(?7)],  (9) 

169  where  X  and  p  are  the  upward  and  downward  transition  rates,  respectively.  It  should  be  noted 

170  that  these  transitions  are  small  and  the  parameters  have  been  carefully  chosen  such  that  the  long- 

171  term  average  of  the  stochastic  model  solutions  recovers  the  dynamics  of  the  original 

172  deterministic  model.  Additional  details  regarding  the  stochastic  version  of  the  skeleton  model  are 

173  provided  in  Thual  et  al.  [2014]. 

174  This  study  continues  to  apply  the  previous  condition  where  the  external  sources  of 

175  cooling  and  moistening  are  set  equal  (i.e.,  sq  =  s°),  which  conserves  the  vertically  integrated 

176  moist  static  energy.  This  balance  condition  was  used  for  both  a  horizontally  homogeneous  and 

177  idealized  warm  pool  state  in  Majda  and  Stechmann  [2009,  201 1]  and  Thual  et  al.  [2014],  where 

178  the  latter  is  shown  by  a  simple  cosine  pattern  in  Figure  la  with  its  maximum  amplitude  shifted 

179  towards  the  longitudes  where  it  is  actually  observed.  Furthermore,  the  zonal  average  of  the 
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180  external  terms  must  be  equal  (i.e.,  sq  =  s°)  in  order  for  a  steady-state  solution  to  exist  without 

181  damping  in  the  above  model  [e.g.,  Majcla  and  Klein ,  2003].  In  the  interest  of  producing  the  best 

182  possible  comparison  with  the  observed  RMM  index,  we  also  perform  several  new  simulations 

183  using  the  National  Oceanic  and  Atmospheric  Administration  (NOAA)  Optimum  Interpolation 

184  (01)  version  2  (v2)  SSTs  [ Reynolds  el  al.,  2002]  in  order  to  provide  a  more  realistic  estimate  of 

185  the  zonal  variability  of  the  external  forcing  terms  used  in  the  model  simulations.  While  the  above 

186  approach  is  far  simpler  than  using  observed  estimates  of  diabatic  heating,  convective  heating, 

187  and  latent  heat  fluxes  to  force  the  skeleton  model  [Ogrosky  and  Stechmann,  2015],  it  is 

188  nevertheless  a  meaningful  assumption  that  zonal  variations  in  moisture  anomalies,  radiative 

189  cooling,  and  convection  arise  primarily  from  the  distribution  of  SSTs,  particularly  for  a  long- 

190  term  climatology  as  considered  in  this  study. 

191  The  annual  average  NOAA  OI  SSTs  v2  (1971-2000)  were  meridionally  averaged  from 

192  10°S-10°N  and  regridded  to  match  the  relatively  coarse  horizontal  resolution  of  the  skeleton 

193  model  (625  km).  A  Poisson  relaxation  was  also  performed  to  interpolate  over  missing  data  (i.e., 

194  land)  and  the  resulting  zonal  profile  was  smoothed  such  that  it  only  retains  contributions  from 

195  the  first  eight  Fourier  modes  (Figure  lb).  Overall,  the  observed  SSTs  closely  match  the  shape  of 

196  the  idealized  warm  pool,  especially  over  the  Maritime  Continent  and  Pacific  Ocean  (120°E- 

197  120°W).  Most  notably,  there  is  a  secondary  maximum  in  the  observed  SSTs  over  the  Atlantic 

198  Ocean  (60°W-0°E)  and  a  slightly  weaker  zonal  gradient  over  the  Indian  Ocean  (60°E-120°E) 

199  than  otherwise  present  in  the  idealized  warm  pool  case.  The  effect  of  these  differences  on  the 

200  MJO-like  waves  in  the  skeleton  model  will  be  further  explored  in  section  3. 

201 

202  2.2.  Defining  an  RMM -like  Index 
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203  As  previously  mentioned,  Stechmann  and  Majda  [2015]  developed  an  MJO  skeleton 

204  index  that  successfully  identified  the  MJO  in  observational  data  by  projecting  reanalysis  onto  the 

205  MJO  modes  obtained  from  the  skeleton  model  when  using  a  uniform  background  state  [i.e.,  sq(x) 

206  =  s9(x)  =s (>\ .  Thus,  rather  than  an  empirical  basis  function  like  the  RMM  indices,  this  index  uses 

207  a  theoretically  based  basis  function  obtained  from  the  model.  This  analysis  was  extended  in 

208  Ogrosky  and  Stechmann  [2015]  to  a  case  with  zonally  varying  forcing  where  sq(x)  ^  s°(x),  which 

209  again  summarizes  the  information  from  all  MJO  wavenumbers  into  a  single  quantity.  While  the 

210  MJO  skeleton  index  has  several  advantages  including  not  requiring  any  time-filtering  nor  the  use 

211  of  EOF  analysis,  its  potential  for  comparison  to  previous  work  is  mostly  limited  due  to  the 

212  relative  newness  of  this  index  combined  with  a  change  in  the  means  of  obtaining  the  basis 

213  function. 

214  In  the  present  work,  we  develop  an  MJO  index  for  the  skeleton  model  that  to  first  order 

215  contains  many  similarities  to  the  RMM  index  of  Wheeler  and  Hendon  [2004].  Although  there  are 

216  many  other  indices  that  are  currently  used  to  identify  the  MJO  including  wavenumber-frequency 

217  filtering  [e.g.,  Wheeler  el  al.,  2000],  lower-tropospheric  wind  [e.g.,  Mcdoney  and  Hartmann , 

218  1998],  upper-troposphere  divergence  and  velocity  potential  [e.g.,  Chen  and  Del  Genio,  2009; 

219  Ventrice  et  al,  2013],  or  time-filtered  outgoing  longwave  radiation  (OLR)  and  cloudiness 

220  datasets  [e.g.,  Hendon  and  Scdby,  1994;  Matthews ,  2000;  Matthews ,  2008;  Kiladis  et  al.,  2014], 

22 1  the  choice  of  the  RMM  index  is  particularly  appropriate  here  given  the  fact  that  the  RMM  index 

222  is  primarily  sensitive  to  wavenumber- 1  modes,  planetary  scale  circulation  anomalies,  and  a 

223  baroclinic  structure  [e.g.,  Straub,  2013]  that  are  likely  to  be  reproduced  by  the  skeleton  model 

224  given  its  focus  on  these  scales. 
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225  The  RMM  is  a  combined  index  of  cloudiness  and  circulation  anomalies  based  on 

226  observed  OLR  and  lower-  and  upper-tropospheric  zonal  winds  taken  from  reanalysis  data  (U850 

227  and  U200,  respectively).  Accordingly,  we  must  first  decide  what  model  variables  are  most 

228  appropriate  to  use  as  observational  proxies  when  creating  an  RMM-like  index  from  the  skeleton 

229  model.  The  following  analysis  uses  the  zonal  wind  from  the  skeleton  model  as  a  direct  substitute 

230  for  the  lower-tropospheric  wind,  whereas  the  convective  (latent)  heating  variable  Ha  is  used  as  a 

23 1  proxy  for  the  observed  OLR.  Stechmann  and  Ogrosky  [2014]  showed  that  variations  in  the  total 

232  diabatic  heating  (i.e.,  Ha  —  s°)  were  directly  proportional  to  OLR  in  the  tropics,  especially  when 

233  averaged  for  long  time  periods  (i.e.,  multiple  years  or  decades)  as  done  in  our  analysis.  As  such, 

234  positive  Ha  anomalies  (enhanced  heating)  can  thus  be  thought  to  have  the  same  physical 

235  interpretation  as  negative  OLR  anomalies  (enhanced  cloudiness  and  deep  convection)  in  our 

236  index. 

237  Although  the  RMM  index  uses  three  variables  (OLR,  U850,  and  U200),  the  choice  of  a 

238  bivariate  index  does  not  preclude  a  meaningful  comparison  of  the  skeleton  model  to 

239  observations.  Straub  [2013]  calculated  time  series  of  the  RMM  principal  components  using  a 

240  method  identical  to  Wheeler  and  Hendon  [2004]  with  additional  experiments  involving  data 

241  withholding.  Straub  [2013]  found  that  the  bivariate  correlations  between  the  full  RMM  index  and 

242  RMM  indices  calculated  without  U850  (i.e.,  OLR  and  U200  only)  or  U200  (i.e.,  OLR  and  U850) 

243  were  0.95  and  0.94  for  the  entire  1979-2010  period,  respectively.  Furthermore,  time  series 

244  constructed  from  a  univariate  analysis  when  considering  only  lower-  or  upper-troposphere  zonal 

245  winds  had  correlation  values  in  excess  of  0.9  with  the  full  RMM.  Consequently,  the  RMM  index 

246  is  largely  dependent  on  the  circulation  component  alone  and  using  a  single-level  for  the  zonal 

247  wind  as  done  here  should  be  sufficient  for  the  detection  and  monitoring  of  planetary  scale  MJOs. 
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248  Following  Wheeler  and  Hendon  [2004],  daily  anomalies  of  the  skeleton  model  variables 

249  were  initially  calculated  by  removing  the  long-term  mean  and  first  three  harmonics  of  the  annual 

250  cycle.  It  should  be  noted  that  while  Wheeler  and  Hendon  [2004]  considered  global  data  between 

25 1  15°S-15°N,  the  skeleton  model  analysis  is  restricted  to  data  at  the  equator  since  the  initial 

252  truncation  only  retains  the  first  meridional  mode  (i.e.,  no  asymmetric  modes).  Wheeler  and 

253  Hendon  [2004]  continue  to  remove  longer-term  variability  by  subtracting  the  component  linearly 

254  related  to  a  measure  of  ENSO  SSTs  and  further  subtract  a  120-day  mean  of  the  previous  120 

255  days  to  remove  any  additional  aspects  of  interannual  or  decadal  variability  and  long-term  trends. 

256  Since  the  skeleton  model  does  not  account  for  ENSO  physics  and  variability  nor  contains  any 

257  time-dependent  SSTs  (recall,  sq  and  s°  are  fixed),  we  choose  to  instead  time  filter  the  daily 

258  anomalies  using  a  20-100  day  Lanczos  filter  with  201  weights  in  order  to  better  isolate  the 

259  intraseasonal  signal  in  the  model  following  Matthews  [2008].  Each  variable  in  the  time-filtered 

260  dataset  is  normalized  by  its  global  variance  to  ensure  equal  contributions  to  the  combined  EOF 

261  analysis,  again  using  the  methods  documented  in  Wheeler  and  Hendon  [2004].  Both  univariate 

262  and  multivariate  EOF  analyses  were  performed  on  the  convective  heating  and  zonal  wind  data 

263  taken  from  the  skeleton  model  simulations.  The  filtered  model  data  were  then  projected  onto  the 

264  leading  two  EOF  modes  to  calculate  the  model  principal  components.  Finally,  the  principal 

265  components  from  each  simulation  and  corresponding  EOF  analysis  were  normalized  by  their 

266  standard  deviation  such  that  a  value  of  unity  represents  anomalies  of  one  standard  deviation  from 

267  the  mean. 

268  We  refer  to  the  RMM-like  index  for  the  model  simulations  as  the  skeleton  multivariate 

269  MJO  index  (SMM)  as  we  only  briefly  consider  those  cases  using  a  univariate  analysis.  Likewise, 

270  the  principal  components  from  the  skeleton  model  data  are  termed  SMM1  and  SMM2  and  are 


13 


27 1  assigned  their  order  based  on  which  of  the  two  leading  EOFs  most  closely  resembles  the 

272  longitudinal  structures  of  the  observed  RMM  index  shown  in  Wheeler  and  Hendon  [2004].  The 

273  values  of  SMM1  and  SMM2  can  therefore  be  visualized  in  a  2-dimensional  phase  space  similar 

274  to  the  RMM  index  with  significant  MJO  activity  defined  as  an  SMM  amplitude  greater  than  or 

275  equal  to  one  with  eastward  propagation  indicated  by  a  counter-clockwise  rotation  of  the  daily 

276  amplitudes  in  the  corresponding  phase  space. 

277 

278  3.  Results 

279  The  following  analysis  presents  a  comparison  of  the  skeleton  model  simulations  with 

280  either  an  idealized  warm  pool  or  the  truncated,  observed  SSTs  using  the  model  set-up  and 

281  framework  established  by  Thual  el  al.  [2014].  The  default  parameters  and  model  constants  for  all 

282  simulations  were  also  taken  from  Thual  et  al.  [2014],  with  the  sole  exception  that  the  zonal 

283  average  of  the  radiative  cooling  term,  se(x),  is  set  equal  to  1.0  K  day"1  and  1.33  K  day'1  for  the 

284  idealized  warm  pool  and  observed  SST  cases,  respectively.  Although  the  skeleton  model  has 

285  recently  been  amended  to  include  a  seasonally  migrating  warm  pool  and  multiple  meridional 

286  modes  [ Thual  et  al.,  2015],  in  addition  to  a  refined  vertical  structure  [ Thual  and Majda,  2015], 

287  we  restrict  our  initial  discussion  to  the  stochastic  and  deterministic  model  simulations  using  the 

288  idealized  and  observed,  annual  average  SSTs. 

289 

290  3.1.  Climatology  and  Variance  of  Simulations  using  Idealized  and  Obserx’ed  Forcing 

291  Time-longitude  series  of  the  anomalous  20-100  day  filtered  skeleton  model  variables 

292  (Ha,  u,  and  q)  are  shown  in  Figure  2  for  a  representative  600  day  period  for  the  simulations 

293  using  the  idealized  warm  pool  and  truncated,  observed  SSTs  (Case-WP  and  Case-SST8N, 
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294  respectively).  In  both  cases,  there  appear  to  be  several  instances  of  eastward  propagating,  MJO- 

295  like  waves  that  travel  at  approximately  5-6  m  s"1  throughout  the  domain.  These  waves  are  evident 

296  throughout  the  entire  period  for  Case-WP  in  the  zonal  wind  (Figure  2b)  with  coherent  eastward 

297  propagating  structures  not  appearing  until  after  day  38200  for  Case-SST8N  (Figure  2e). 

298  Although  specific  MJO  events  should  not  be  directly  compared  between  the  two  simulations, 

299  Case-SST8N  displays  greater  variance  in  the  convective  heating  over  the  western  Indian  Ocean 

300  (i.e.,  60-90°E)  and  a  stronger  mean  convective  heating  over  the  tropical  Atlantic  (60°W-0°E) 

301  than  identified  for  Case-WP  (cf.  Figures  2a  and  2d).  These  areas  of  enhanced  heating  agree 

302  exceedingly  well  with  Ogrosky  and  Stechmann  [2015],  who  calculated  observed  Ha  values  from 

303  precipitation  data  using  the  Global  Precipitation  Climatology  Project  (GPCP)  [Huffman  et  al., 

304  2012]  and  similar  model  quantities  from  simulations  forced  with  observed  atmospheric  OLR  and 

305  latent  heat  fluxes  derived  from  reanalysis  data  [Figure  13  in  Ogrosky  and  Stechmann ,  2015]. 

306  Power  spectra  of  the  reconstructed  convective  heating  and  zonal  wind  as  a  function  of  the 

307  zonal  wavenumber  (k)  and  frequency  (co)  are  shown  in  Figure  3  for  a  34  year  period  from  each 

308  model  simulation.  Both  Case-WP  and  Case-SST8N  demonstrate  significant  power  in  the  MJO 

309  range  (i.e.,  1  <  k  <  5  and  1/90  <a><  1/30  cpd)  with  the  peak  power  stretching  along  a  band  with 

310  doo/dk  ~  0,  an  auspicious  and  recurring  feature  of  the  skeleton  model  [Thucd  et  cd.,  2014].  The 

311  peak  power  in  each  case  is  also  near  the  dispersion  curves  for  the  linear  solutions  from  the 

312  skeleton  model  (the  black  circles  in  Figure  3).  As  in  Wcdiser  et  cd.  [2009],  the  breadth  of  the  peak 

313  in  frequency  and  wavenumber  for  the  convective  heating  (i.e.,  a  proxy  for  precipitation)  relative 

314  to  zonal  wind  is  captured  by  the  skeleton  model,  though  the  convective  heating  power  in  the 

315  MJO  band  is  nearly  constant  with  increasing  wavenumber  unlike  in  observations  where  the 

316  precipitation  power  is  largely  confined  to  1  <  k  <  3.  Generally,  the  convective  heating  and  zonal 
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317  wind  for  Case-SST8N  shows  increased  power  near  the  dispersion  curves  of  both  the  MJO  and 

318  moist  Rossby  modes,  likely  due  to  the  slightly  higher  values  of  s9  and  se in  this  region  (67.5°E) 

319  and  zonal  mean  when  compared  to  the  scenario  with  the  idealized  warm  pool.  The  peak  power 

320  for  the  stochastic  model  simulations  also  occurs  at  slightly  lower  frequencies  than  the  MJO- 

321  mode  dispersion  curves,  similar  to  Thual  et  al.  [2014].  In  addition,  there  are  secondary  peaks  in 

322  Ha  and  u  at  double  and  triple  the  MJO-mode  frequencies,  both  of  which  are  likely  attributed  to 

323  the  nonlinear  qa  term  in  the  skeleton  model. 

324  Figure  4  shows  the  leading  two  structures  from  a  univariate  EOF  analysis  of  the 

325  convective  heating  and  zonal  wind  for  Case-WP.  While  the  variance  in  convective  heating,  Ha, 

326  is  centered  over  the  longitudes  with  the  maximum  forcing  for  the  idealized  warm  pool  case 

327  (Figure  la),  the  EOF  structures  for  the  univariate  analysis  show  significant  noise  due  to  the 

328  presence  of  nonlinear  terms  and  inclusion  of  stochastic  effects  (Figures  4a  and  4b),  even  after 

329  having  already  applied  a  20-100  day  filter.  Furthermore,  EOF1  and  EOF2  only  each  explain 

330  3.4%  of  the  observed  variance  despite  using  a  sufficiently  long  record.  Although  the  skeleton 

33 1  model  is  already  computed  on  a  relatively  coarse  grid  (625  km  or  5.625°)  the  stochastic 

332  processes  in  the  model  occur  at  scales  too  small  to  identify  recurring  modes  with  the  simple 

333  global  EOF  analysis  without  additional  smoothing  first  being  applied  to  the  model  data.  In  fact, 

334  each  of  the  first  eight  modes  explained  approximately  the  same  amount  of  variance  over  the 

335  entire  340  years  of  simulation  data,  indicating  that  the  first  two  EOFs  are  not  unique  (not  shown). 

336  Meanwhile,  the  leading  structures  from  the  univariate  zonal  wind  analysis  for  Case-WP  show  a 

337  highly  symmetric  wavenumber- 1  pattern  centered  about  the  maximum  of  the  idealized  warm 

338  pool  with  31.3%  of  the  variance  explained  by  the  leading  two  modes  (Figures  4c  and  4d).  The 

339  third  and  fourth  EOF  structures  for  Case-WP  identify  a  wavenumber-2  circulation  anomaly 
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340  centered  about  the  warm  pool  maximum  that  explain  an  additional  16.2%  of  the  combined 

341  variance  (not  shown). 

342  The  first  two  modes  from  the  multivariate  EOF  analysis  for  Case-WP  are  shown  in 

343  Figures  5a  and  5b.  The  order  of  the  leading  EOFs  has  been  switched  to  better  facilitate 

344  comparisons  with  the  observed  EOF  patterns  from  Wheeler  and  Hendon  [2004],  which  are 

345  plotted  for  reference  in  Figures  5c  and  5d.  The  EOF  structures  for  the  convective  heating 

346  anomalies  improve  significantly  when  constrained  by  the  variations  in  zonal  wind  for  the 

347  combined  analysis  with  a  clear  zonal  wavenumber- 1  pattern  of  enhanced  heating  and  cooling  (cf. 

348  Figures  4a-4b  and  Figures  5a-5b).  Albeit,  the  previous  improvement  comes  at  the  expense  of  a 

349  reduction  in  the  total  explained  variance,  which  is  only  16.7%  for  the  combined  EOF  analysis. 

350  The  heating  EOFs  have  a  maximum  correlation  coefficient  of  0.70  and  0.68  with  the  observed 

35 1  Wheeler  and  Hendon  [2004]  OFR  EOFs  (Figures  5c  and  5d)  when  calculated  at  all  longitudes  in 

352  order  to  account  for  any  longitudinal  displacement  of  the  convective  features  in  the  model  (Table 

353  1).  The  zonal  wind  structures  in  the  multivariate  analysis  are  nearly  identical  to  the  univariate 

354  patterns,  suggesting  that  the  combined  index  is  far  more  sensitive  to  the  circulation  component 

355  alone  as  is  the  case  in  observations  [e.g.,  Straub ,  2013].  The  multivariate  circulation  modes 

356  generally  mimic  the  structures  in  the  observed  Wheeler  and  Hendon  [2004]  EOFs  with  a 

357  maximum  correlation  of  0.83  and  0.72  (Table  1),  though  the  observed  EOF  patterns  contain 

358  more  asymmetry  and  explain  a  higher  percent  of  the  total  variance  (25.0%)  than  Case-WP  from 

359  the  skeleton  model.  Regardless,  it  is  impressive  the  skeleton  model  is  able  to  reproduce  the 

360  salient  features  of  the  observed  RMM  index  considering  the  many  simplifications  in  the  model 

361  and  use  of  an  idealized  warm  pool  in  lieu  of  observed  SSTs. 
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362  The  equivalent  multivariate  EOF  modes  for  the  skeleton  model  simulation  using  the 

363  truncated,  observed  SSTs  (Case-SST8N)  are  shown  in  Figures  6a  and  6b.  Compared  to  Case- 

364  WP,  the  addition  of  the  observed  sq  and  .^  increases  the  asymmetry  in  the  convective  heating  and 

365  zonal  wind  structures  and  better  matches  the  observed  Wheeler  and  Hendon  [2004]  EOFs 

366  (Figures  6c  and  6d).  As  such,  the  maximum  correlation  coefficient  between  the  model  and 

367  observed  EOFs  increases  from  0.70  and  0.68  (Case-WP)  to  0.89  and  0.80  (Case-SST8N)  for  the 

368  convective  heating  variable  (Table  1).  Fikewise,  the  correlations  for  the  circulation  anomalies 

369  increase  from  0.83  and  0.72  to  an  incredible  0.98  and  0.94  for  Case-WP  and  Case-SST8N, 

370  respectively  (Table  1)  As  in  Case-WP,  the  model  EOF  structures  are  in  quadrature  with  the 

37 1  maximum  cross-correlation  (i.e.,  between  the  corresponding  principal  components)  for  Case- 

372  SST8N  at  9  days  (Figure  7),  consistent  with  observations  and  suggesting  an  overall  MJO 

373  oscillation  frequency  of  ~40  days. 

374  The  peak  positive  values  in  zonal  wind  for  the  Case-SST8N  multivariate  EOFs  are 

375  located  about  30-40°  east  of  where  the  maxima  occur  in  observations  (Figure  6).  The  combined 

376  EOFs  also  explain  15.7%  of  the  total  variance,  which  is  slightly  lower  than  for  Case-WP. 

377  Nevertheless,  the  multivariate  EOFs  bear  a  striking  resemblance  to  the  observed  EOF  patterns 

378  from  the  RMM  index,  with  maximum  correlations  of  0.98  and  0.94  for  the  circulation 

379  component  and  0.89  and  0.80  for  the  convection  anomalies  (Table  1).  As  such,  we  therefore 

380  choose  to  also  project  the  model  data  directly  onto  the  Wheeler  and  Hendon  [2004]  EOFs  when 

381  deriving  various  forms  of  the  model  principal  components  (i.e.,  hybrid  SMM1  and  SMM2  values 

382  that  depend  on  both  the  model  and  observations)  in  order  to  directly  evaluate  how  well  the 

383  skeleton  model  is  able  to  reproduce  observed  MJO  event  statistics  using  the  RMM  framework. 

384  The  use  of  a  reference  set  of  EOFs  for  calculating  the  model  principal  components  also  provides 
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385  a  valuable  standard  and  better  facilitates  comparison  of  MJO  climatology  statistics  from  different 

386  model  simulations  using  varying  parameters  or  forcing  data. 

387 

388  3.2.  Identifying  MJO  Events  in  the  Skeleton  Model 

389  Having  already  developed  an  RMM-like  index  that  is  similar  to  observations,  we  now 

390  turn  our  attention  to  applications  of  the  SMM  index,  including  the  identification  of  MJO  events 

391  in  the  skeleton  model.  Similarly,  we  also  document  the  sensitivities  of  using  the  model-derived 

392  univariate  and  multivariate  EOFs,  in  addition  to  the  Wheeler  and  Hendon  [2004]  observed 

393  multivariate  EOFs  for  the  resulting  MJO  event  climatology. 

394  This  study  uses  the  common  practice  of  one  standard  deviation  (i.e.,  a  combined 

395  amplitude  greater  than  or  equal  to  unity)  as  a  lower  threshold  for  defining  the  presence  of  a 

396  strong  MJO.  In  addition,  we  employ  the  technique  used  in  Stachnik  et  al.  [2015],  where  the  daily 

397  values  of  RMM1  and  RMM2  must  contain  nearly  continual  counter-clockwise  rotation  (i.e., 

398  eastward  propagation)  in  the  2-dimensional  phase  space  of  Wheeler  and  Hendon  [2004].  We  also 

399  require  that  a  simulated  MJO  “event”  pass  through  at  least  four  phases  of  the  SMM  index. 

400  Two  examples  of  simulated  MJO-like  waves  are  shown  in  Figure  8  for  Case-WP.  For 

401  each  event,  the  SMM  index  was  calculated  by  projecting  the  skeleton  model  data  onto  the 

402  multivariate,  observed  Wheeler  and  Hendon  [2004]  EOFs.  Using  the  above  definition,  a  primary 

403  MJO  event  initiates  on  day  30564  and  day  32665  for  the  periods  in  Figures  8a  and  8b  and  are 

404  noted  with  a  gray  circle  for  the  first  days  when  the  combined  SMM  amplitude  grows  above  1.0 

405  along  the  far  eastern  edge  of  the  RMM  phase  8  domain  and  western  region  of  RMM  phase  6, 

406  respectively.  The  MJO  event  in  Figure  8a  begins  to  propagate  eastward  through  two  phases 

407  (RMM  phase  8  and  phase  1)  before  being  identified  as  a  continuing  MJO  event  (in  contrast  to 
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408  primary  events  that  develop  from  inactive  conditions)  provided  that  it  propagate  through  four 

409  additional  phases  as  is  highlighted  in  blue.  These  conditions  are  met  again  as  the  MJO  moves 

410  into  the  RMM  phase  3  and  phase  4  domains,  each  counting  as  a  separate  continuing  event  with 

411  the  subsequent  paths  colored  in  green  and  orange,  respectively.  For  simplicity,  only  the  first  day 

412  of  the  continuing  events  are  shown  in  Figure  8b,  with  a  special  circumnavigating  event  or  wave 

413  train  MJO  denoted  in  gold  for  RMM  phase  6  once  the  MJO-like  wave  has  successfully 

414  propagated  around  all  longitudes.  As  such,  the  initiation  of  each  primary  MJO  may  contain 

415  multiple  continuing  events,  especially  for  very  long-lasting  MJOs. 

416  Termination  events  occur  when  the  SMM  amplitude  decays  below  unity  for  an  extended 

417  period  of  time  or  demonstrates  significant  westward  propagation.  Examples  of  MJO  termination 

418  events  are  shown  in  RMM  phase  7  and  phase  1  for  Figures  8a  and  8b,  respectively.  Overall,  the 

419  event  histories  are  much  smoother  than  in  observations  [e.g.,  Figure  1  in  Stachnik  el  al.,  2015], 

420  even  when  projecting  the  model  data  directly  onto  the  Wheeler  and  Hendon  [2004]  EOFs.  This 

42 1  artifact  is  partially  by  design,  as  the  skeleton  model  was  specifically  developed  for  intraseasonal 

422  and  planetary  oscillations.  The  use  of  a  band-pass  filter  on  the  model  data,  while  able  to  better 

423  isolate  the  intraseasonal  signal,  lacks  the  same  day-to-day  noise  and  high  frequency  variability 

424  originating  from  convectively  coupled  equatorial  waves  and  other  disturbances  that  are  capable 

425  of  projecting  upon  the  RMM  index  in  nature  [e.g.,  Roundy  et  al.,  2009].  Additional  details  and 

426  sensitivities  of  the  search  algorithm  are  described  in  Stachnik  et  al.  [2015]. 

427  Table  2  lists  the  number  of  MJO  events  from  34  years  of  observations  (1979-2012)  using 

428  the  historical  RMM  index  and  the  mean  and  95%  confidence  intervals  for  an  equivalent  34  year 

429  period  derived  from  680  years  of  skeleton  model  data  using  an  idealized  warm  pool  (Case-WP). 

430  Overall,  the  number  of  simulated  MJO  events  is  in  fairly  good  agreement  with  observations 
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43 1  when  one  considers  the  highly  idealized  form  of  the  model  and  the  lack  of  model  tuning  relative 

432  to  this  metric.  In  detail,  for  all  methods  of  calculating  the  SMM  index  (i.e.,  univariate  and 

433  multivariate  model  EOFs  in  addition  to  the  multivariate  observed  EOFs)  the  number  of  primary 

434  MJOs  is  less  than  observations  with  a  range  89.3-97.7  events  for  the  skeleton  model  compared  to 

435  154  events  in  observations.  The  uncertainty  estimates  for  primary  MJO  events  are  fairly  small, 

436  with  a  limited  range  of  3. 6-4. 8%  relative  to  the  mean  number  of  events.  With  the  exception  of 

437  the  univariate  zonal  wind  EOF  analysis,  the  mean  number  of  continuing  events  exceeds 

438  observations,  though  the  observed  number  falls  within  the  95%  confidence  intervals  in  each  case. 

439  Similarly,  all  of  the  Case-WP  scenarios  contain  approximately  40-60%  more  circumnavigating 

440  events  than  actually  occur.  Finally,  the  number  of  terminating  events  is  identical  to  the  primary 

441  MJOs,  as  all  events  that  initiate  must  eventually  come  to  an  end  in  the  skeleton  model. 

442  The  decrease  in  the  number  of  separate  MJO  events  in  the  skeleton  model  while 

443  maintaining  a  nearly  equal  or  slightly  greater  number  of  continuing  and  circumnavigating  events 

444  implies  that  the  average  MJO  event  length  (i.e.,  the  number  of  days  between  event  initiation  and 

445  termination)  should  be  higher  in  the  model  than  observations.  This  presumption  is  correct,  with  a 

446  mean  model  event  length  of  about  40  days  for  all  realizations  from  Case-WP  in  contrast  to  an 

447  approximate  35  day  event  length  in  observations  (Table  2).  Also,  it  is  interesting  to  note  that 

448  there  is  little  overall  sensitivity  to  the  choice  of  EOFs  when  creating  the  model  event 

449  climatology.  Nevertheless,  using  the  observed  Wheeler  and  Hendon  [2004]  EOFs  resulted  in  the 

450  most  favorable  comparison  for  primary  MJOs  with  those  seen  in  nature. 

45 1  Similar  results  for  the  number  of  MJO  events  and  average  MJO  event  length  are  provided 

452  in  Table  3  for  the  case  using  the  truncated,  observed  SSTs  (Case-SST8N).  Once  again,  the 

453  number  of  primary  MJOs  is  less  than  in  observations,  with  only  about  64%  the  number  of 
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454  observed  MJOs  regardless  of  the  choice  in  EOFs.  Here,  the  percent  increase  in  the  number  of 

455  continuing  and  circumnavigating  events  is  higher  than  for  Case-WP,  with  443.5-477.9 

456  continuing  events  and  33.4-36.6  circumnavigating  MJOs  in  the  skeleton  model  compared  to  330 

457  and  15  MJO  events  in  observations,  respectively.  As  a  result,  the  average  MJO  event  length  for 

458  Case-SST8N  is  higher  than  the  idealized  warm  pool  simulation,  with  an  average  length  of 

459  approximately  44  days. 

460  The  increase  in  the  average  MJO  event  length  for  the  skeleton  model  simulations  is 

46 1  somewhat  deceptive,  however,  as  further  investigation  revealed  that  the  mean  event  length  is 

462  primarily  skewed  by  a  smaller  number  of  simulated  events  that  have  very  long  length.  Figure  9 

463  shows  the  cumulative  probability  distribution  function  of  the  observed  (black)  and  modeled 

464  (blue)  MJO  event  lengths  for  a  representative  34  year  simulation  period  using  the  truncated, 

465  observed  SSTs.  Approximately  two-thirds  of  the  observed  events  fall  within  20-100  days  with  a 

466  slightly  lower  (-60%)  percent  of  simulated  MJOs  from  Case-SST8N  in  the  same  range  (Figure 

467  9).  Furthermore,  the  shortest  40-50%  of  simulated  MJO  events  in  the  stochastic  version  of  the 

468  skeleton  model  match  observations. 

469  Moreover,  the  introduction  of  stochasticity  in  the  model  is  critical  for  generating  a  more 

470  realistic  MJO  event  climatology.  Additional  simulations  using  the  deterministic  version  of  the 

47 1  model  with  the  idealized  warm  pool  or  observed  SSTs  (Case-WPdet  and  Case-SST8Ndet, 

472  respectively)  showed  about  a  10%  reduction  in  the  number  of  primary  MJOs  compared  to  their 

473  corresponding  stochastic  simulation.  The  deterministic  solutions  also  had  an  approximate  130% 

474  and  175%  increase  in  the  number  of  continuing  events  for  Case-WPdet  and  Case-SST8Ndet, 

475  respectively  (Table  4).  The  resulting  average  MJO  event  lengths  for  the  deterministic  solutions 

476  far  exceeds  observations,  with  Figure  9  indicating  a  significant  positive  bias  for  all  MJO  event 
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All  lengths  for  a  representative  34  year  period  of  Case-SST8Ndet  (red).  In  fact,  nearly  10%  of  the 

478  simulated  MJOs  last  more  than  150  days  in  this  case. 

479 

480  3.3.  Geographic  Distribution  ofMJO  Events  in  the  Skeleton  Model 

48 1  Thus  far,  the  use  of  real  data  SSTs  in  the  skeleton  model  has  been  shown  to  produce  more 

482  realistic  multivariate  EOF  structures  with  higher  correlation  values  that  better  match  the 

483  observed  Wheeler  and  Hendon  [2004]  RMM  EOFs.  Furthermore,  those  model  simulations  using 

484  observed  forcing  [Ogrosky  and  Stechmann,  2015]  produce  a  more  realistic  geographic 

485  distribution  of  convective  heating  with  an  increase  in  Ha  variance  over  the  locations  where  it  is 

486  actually  observed.  Therefore,  it  becomes  reasonable  to  wonder  how  zonal  variations  in  the 

487  external  forcing  terms  might  alter  the  geographic  distribution  of  MJO  initiation  and  termination 

488  events  in  the  skeleton  model.  Prior  to  applying  this  test  for  natural  climate  variability  such  as  the 

489  modulated  forcing  due  to  ENSO  or  Indian  Ocean  dipole  (IOD)  SSTs,  it  is  first  necessary  to 

490  determine  whether  the  skeleton  model  can  reproduce  the  observed  geographic  variability  of  MJO 

491  initiation  and  termination  events  using  annual  average  SSTs. 

492  Figure  10  shows  the  relative  frequency  of  occurrence  for  primary,  continuing,  and 

493  terminating  MJO  events  subset  by  RMM  phase  for  34  years  of  observations  and  similar  periods 

494  of  model  data  for  Case-WP,  Case-SST8N,  and  Case-SST8Ndet  using  the  observed  Wheeler  and 

495  Hendon  [2004]  EOFs  when  calculating  the  model  climatology.  While  the  figures  are  similar  to 

496  the  2-dimensional  phase  space  diagrams  from  Wheeler  and  Hendon  [2004],  the  plots  indicate  the 

497  percent  likelihood  of  each  MJO  event  type  visualized  in  a  radial-frequency  domain.  On  average, 

498  observed  MJO  events  are  most  likely  to  initiate  over  the  Indian  Ocean  (RMM  phase  2)  and 

499  terminate  over  the  Western  Pacific  or  western  hemisphere  (RMM  phase  6  and  phase  8),  as 
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500  evidenced  in  Figure  10a.  Figure  10b  shows  an  analogous  distribution  for  a  single  34  year  period 

501  for  Case-WP.  There  is  some  suggestion  that  the  simulated  MJOs  prefer  to  initiate  in  the  RMM 

502  phase  1  and  phase  2  domains,  with  local  termination  maxima  in  RMM  phase  8  and  phase  5  (i.e., 

503  western  hemisphere  and  exiting  the  Maritime  Continent,  respectively).  Comparisons  for  Case- 

504  SST8N  (Figure  10c)  are  in  less  agreement,  however,  with  a  large  termination  peak  located  within 

505  the  RMM  phase  3  domain  that  does  not  match  observations.  The  deterministic  model  simulation 

506  for  this  case  (Case-SST8Ndet)  is  also  different  from  observations  with  the  simulated  MJOs 

507  becoming  too  regular  and  preferentially  beginning  and  terminating  in  RMM  phase  3  and  phase  2, 

508  respectively  (Figure  lOd). 

509  To  better  determine  whether  or  not  the  skeleton  model  is  capable  of  reproducing  the 

510  observed  MJO  event  climatology  as  a  function  of  the  large-scale  forcing,  we  extend  the  previous 

5 1 1  analysis  to  680  years  of  model  data,  i.e.,  20  independent  samples  of  34  year  periods  similar  in 

512  length  to  the  observations.  Figure  1 1  displays  the  mean  event  frequencies  and  95%  confidence 

513  intervals  as  a  function  of  RMM  phase  for  the  observations  and  five  different  model  simulations. 

514  In  all  panels,  the  uncertainty  is  calculated  using  a  binomial  proportion  confidence  interval  that  is 

5 1 5  dependent  on  the  pass  and  fail  rates  of  the  N  events  being  assigned  to  a  particular  RMM  phase 

516  throughout  the  entire  record.  As  such,  the  results  are  only  statistically  significant  if  the 

517  confidence  intervals  do  not  overlap  the  reference  line  of  12.5%,  i.e.,  the  equal  likelihood  that  an 

518  MJO  event  will  occur  within  any  of  the  eight  RMM  phases. 

519  Unfortunately,  most  of  the  frequencies  in  Figure  11a  cannot  be  treated  as  statistically 

520  significant  or  different  from  an  equal  chance  due  to  the  relatively  large  confidence  intervals  and 

52 1  shorter  historical  record  for  observed  MJOs.  The  skeleton  model,  however,  is  not  constrained  by 

522  the  lack  of  samples  and  can  be  run  for  extended  periods.  As  a  check,  we  perform  the  same 
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523  analysis  on  a  special  long-term  stochastic  simulation  using  a  homogenous  background  (Case- 

524  HG)  where  the  principal  components  are  calculated  using  the  model’s  own  multivariate  EOFs, 

525  while  the  other  simulations  use  the  observed  Wheeler  and  Hendon  [2004]  EOFs.  This  situation 

526  should  be  most  likely  to  produce  equal  chance  results  for  all  RMM  phases  given  the  absence  of 

527  geographically  varying  forcing  or  data  projection  upon  the  observed  EOFs  and  is  successfully 

528  recovered  in  Figure  1  Id.  Only  a  few  of  the  mean  frequencies  for  Case-WP  or  Case-SST8N 

529  (Figures  1  lb  and  1  lc)  demonstrate  any  statistical  significance  from  an  equal  probability  and  the 

530  peaks  in  the  MJO  event  distributions  do  not  necessarily  match  observations.  It  remains  possible, 

53 1  however,  that  the  stochastic  effects  in  the  model  act  to  dilute  the  geographic  dependencies  and 

532  impacts  of  applying  zonally  varying  forcing.  The  long-term  climatology  for  Case-WPdet  (Figure 

533  1  le)  does  show  some  features  in  common  with  the  observed  locations  of  MJO  event  initiation 

534  that  are  statistically  significant  yet  demonstrate  little  skill  for  the  corresponding  climatology  of 

535  termination  events.  The  frequency  distribution  for  Case-SST8Ndet  (Figure  1  If)  produces  the 

536  largest  differences  from  an  equal  chance,  though  those  phases  with  distinguishable  frequencies 

537  do  not  match  observations.  As  such,  more  work  is  needed  to  help  separate  the  effects  of 

538  nonlinearity,  stochasticity,  and  zonally-varying  forcing  upon  where  MJO  events  begin  and  decay 

539  in  the  skeleton  model.  Similarly,  the  model  may  need  to  account  for  other  factors  such  as 

540  interactions  with  the  extratropical  circulation  before  making  any  firm  conclusions  regarding  its 

541  potential  future  use  in  studying  natural  climate  variability  and  change. 

542 

543  3.4.  Composites  and  Precursors  for  MJO  Initiation  and  Termination 

544  Figure  12  shows  the  composite  RMM1-RMM2  indices  associated  with  primary  MJO 

545  event  initiation  for  all  phases  in  observations  and  the  corresponding  SMM1-SMM2  histories 
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546  from  a  single  34  year  period  of  skeleton  model  data  for  Case-WP,  Case-SST8N,  and  Case- 

547  SST8Ndet.  The  observed  Wheeler  and  Hendon  [2004]  EOFs  were  used  to  calculate  the  principal 

548  components  in  each  model  case  and  the  day  0  (i.e.,  first  day  with  combined  amplitude  greater 

549  than  1.0)  is  noted  with  a  circumscribed  X  in  all  panels. 

550  In  the  observations  and  stochastic  simulations,  MJO  initiation  occurs  rather  quickly,  with 

551  a  rapid  increase  in  the  combined  amplitudes  for  the  three  days  prior  to  the  initiation  event.  The 

552  observational  composite  shows  continued  growth  for  an  additional  2-3  days  with  most  phases 

553  achieving  a  near  constant  value  with  only  weak  growth  of  the  combined  amplitude  through  the 

554  following  12  days  (Figure  12a).  Meanwhile,  the  stochastic  simulations  share  the  same  initial 

555  growth  (i.e.,  rapid  intensification  for  3  days  prior  to  initiation)  though  continue  to  intensify 

556  beyond  observed  values  and  reach  a  maximum  SMM  amplitude  5-6  days  after  the  initiation  of  an 

557  event  (Figures  12b  and  12c).  Overall,  the  average  root-mean- square  error  (RMSE)  for  all  eight 

558  phases  is  0.27  and  0.23  amplitude  units  for  Case-WP  and  Case-SST8N,  respectively.  The 

559  initiation  composite  for  the  deterministic  simulation  follows  the  same  trend  as  the  observations 

560  and  stochastic  runs  although  the  SMM  amplitude  growth  rates  are  much  slower  before  and  after 

561  MJO  initiation  (Figure  12d).  Consequently,  the  smaller  changes  in  SMM  amplitude  for  Case- 

562  SST8Ndet  result  in  the  largest  values  of  RMSE  (0.29  amplitude  units)  and  do  not  match  the 

563  observations  as  well  at  the  stochastic  simulations. 

564  Analogous  RMM  and  SMM  composites  are  shown  before  and  after  the  start  of  an  MJO 

565  termination  event  in  Figure  13.  Again,  2-3  days  seems  to  be  the  relevant  time  scales  for  the 

566  decay  of  an  observed  MJO  event  as  the  combined  RMM  amplitudes  maintain  nearly  constant 

567  values  in  the  composite  until  approximately  3  days  prior  to  the  termination  markers  that  indicate 

568  the  last  day  where  the  combined  amplitude  remains  above  unity  (Figure  13a).  Only  RMM  phase 
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569  2  termination  events  experience  gradual  decay,  with  the  RMM  amplitude  weakening  for  about  7 

570  days  prior  to  termination  [ Stachnik  et  al.,  2015].  The  composites  for  the  stochastic  model 

571  simulations  attain  stronger  values  of  SMM  amplitude  5-6  days  prior  to  the  eventual  termination 

572  event  than  seen  in  observations  (Figures  13b  and  13c).  The  average  RMSE  values  are  0.27  and 

573  0.26  for  Case-WP  and  Case-SST8N,  respectively.  As  with  initiation,  the  changes  in  SMM 

574  amplitude  for  Case-SST8Ndet  occur  more  slowly  are  much  smaller  than  actually  observed 

575  (Figure  13d).  However,  the  RMSE  for  the  deterministic  run  is  slightly  smaller  (0.24  amplitude 

576  units)  given  the  pronounced  pulsing  in  SMM  amplitude  prior  to  termination  for  the  stochastic 

577  simulations  during  this  period. 

578  In  order  to  identify  possible  mechanisms  for  MJO  termination,  Figure  14  shows  the  lag- 

579  longitude  composites  of  the  anomalous  convective  heating,  moisture,  and  zonal  wind  at  the 

580  equator  from  340  years  of  MJO  events  that  either  continue  or  terminate  in  the  RMM  phase  6 

581  region  for  Case-SST8N  when  using  the  observed  Wheeler  and  Hendon  [2004]  EOFs.  We  focus 

582  on  this  domain  given  the  fact  that  observed  MJOs  often  terminate  in  region  (e.g.,  Figure  10a).  In 

583  addition,  the  RMM6  domain  allows  for  a  longer  analysis  of  the  upstream  conditions  prior  to  the 

584  termination  event  without  the  complication  of  circumnavigating  the  western  hemisphere.  All  of 

585  the  model  variables  have  been  time  filtered  to  isolate  the  20-100  day  intraseasonal  variability. 

586  Differences  between  the  continuing  and  terminating  events  are  shown  in  the  bottom  panel  of 

587  each  column  along  with  stippling  that  indicates  the  presence  of  a  statistically  significant  different 

588  mean  associated  with  continuing  and  terminating  MJOs. 

589  To  begin,  the  differences  in  the  large-scale  environmental  conditions  at  positive  lags  (i.e., 

590  after  the  termination  of  an  MJO  event)  are  not  nearly  as  pronounced  as  they  are  in  observations 

591  [cf.  Figure  4  in  Stachnik  et  al.,  2015].  For  example,  the  zonal  wind  maintains  a  wavenumber-1 
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592  type  circulation  anomaly  after  day  0  for  the  terminating  cases  in  the  skeleton  model  (Figure  14h) 

593  with  the  main  difference  being  a  weaker  amplitude  than  identified  during  positive  lags  for  the 

594  continuing  events  (Figure  14g).  Conversely,  the  wavenumber- 1  circulation  anomaly  decays 

595  rapidly  when  compositing  low-level  zonal  winds  from  reanalysis  data  for  terminating  MJO 

596  events  [Figure  4  in  Stachnik  et  al,  2015].  In  addition,  there  is  no  single  variable  that  seems  to 

597  best  predict  whether  an  MJO  event  will  continue  or  terminate  once  entering  the  RMM  phase  6 

598  domain.  Unlike  in  observations,  significant  differences  between  continuing  and  terminating 

599  events  occur  in  both  the  active  and  inactive  phases  of  the  MJO  for  the  Ha,  q,  and  u  anomalies. 

600  These  difference  appear  to  manifest  within  MJO  wave  trains  and  are  identifiable  more  than  30 

601  days  prior  to  the  eventual  continuing  or  terminating  event.  The  above  result  is  intriguing, 

602  considering  that  the  composite  SMM  indices  in  Figure  13  only  suggest  a  weakening  of  the 

603  SMM1-SMM2  amplitudes  in  the  5-6  days  immediately  prior  to  event  termination.  As  such,  it  is 

604  not  immediately  clear  what  sort  of  physical  precursor  conditions  may  be  affecting  MJO 

605  termination  in  the  skeleton  model. 

606  Figure  15  expands  the  above  analysis  to  include  composites  and  differences  in 

607  convergence  and  potential  temperature  for  the  340  year  period  of  continuing  and  terminating 

608  MJOs.  The  rightmost  column  of  Figure  15  shows  the  potential  temperature  composites  for  a 

609  shorter,  single  34  year  period  of  skeleton  model  data.  Once  more,  the  differences  in  convergence 

610  and  potential  temperature  at  positive  lags  for  termination  events  is  more  closely  related  to 

611  amplitude  differences  rather  than  a  complete  shutdown  of  the  large-scale  pattern.  However,  it  is 

612  interesting  to  note  that  there  are  a  number  of  slow  moving  (~1  m  s"1),  eastward  propagating 

613  waves  embedded  within  the  composite  difference  panels  for  the  convergence  and  potential 

614  temperature  in  Figures  14c,  14f,  15c,  and  15f.  These  anomalies  occur  at  relatively  high 
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615  wavenumber  ( k  ~  12-13)  are  even  more  pronounced  in  the  termination  composites  for  the  shorter 

616  34  year  period  (cf.  Figures  15g  and  15h).  Majda  and  Stechmann  [201 1]  also  noted  potential 

617  interactions  of  dry  Kelvin  and  Rossby  modes  during  MJO  termination  for  the  skeleton  model, 

618  though  the  waves  in  our  analysis  travel  much  more  slowly  than  these  theoretical  modes.  Albeit 

619  highly  speculative,  the  previous  results  suggest  that  continuing  MJO  events  may  suppress  the 

620  development  of  eastward  propagating  waves.  Alternatively,  the  emergence  of  dry  Kelvin  waves 

62 1  in  Majda  and  Stechmann  [2011]  or  other  higher  wavenumber  modes  embedded  within  the  large- 

622  scale  flow  as  seen  here  may  interfere  with  the  MJO  dynamics  and  help  lead  to  its  future  demise. 

623  Finally,  Figure  16  shows  the  distribution  of  zonal  wind  and  moisture  anomalies  for 

624  individual  RMM  phase  6  MJOs  averaged  over  the  local  domain  (120-160°E)  at  10  days  prior  to 

625  the  continuing  or  terminating  event  for  Case-SST8N.  As  is  the  case  with  observations  in 

626  Stachnik  et  cd.  [2015],  there  is  appreciable  overlap  in  the  anomaly  distributions  for  individual 

627  continuing  and  terminating  events  despite  containing  a  statistically  significant  different  mean  for 

628  the  time  and  longitudes  displayed  (e.g.,  Figure  14i).  For  example,  individual  continuing  MJO 

629  events  may  contain  cooling  and  drying  whereas  terminating  events  may  experience  moistening 

630  and  warming  (Figure  16b).  Subsequently,  there  appears  to  be  no  one  recurring  or  unique  set-up 

63 1  for  MJO  termination  even  in  the  relatively  simple  and  low-order  stochastic  skeleton  model. 

632 

633  4.  Summary  and  Discussion 

634  This  study  developed  an  RMM-like  index  for  the  deterministic  and  stochastic  version  of 

635  the  MJO  skeleton  model  [Majda  and  Stechmann,  2009;  201 1]  using  an  approach  similar  to  that 

636  outlined  in  Wheeler  and  Hendon  [2004].  Univariate  and  multivariate  EOF  analyses  were 

637  performed  using  convective  heating  and  zonal  wind  data  taken  from  the  skeleton  model 
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638  simulations  forced  with  an  idealized  warm  pool  and  observed  SSTs.  The  simulations  with  the 

639  observed  SSTs  (Case-SST8N)  demonstrated  a  number  of  small  improvements  over  the  idealized 

640  warm  pool  (Case-WP)  and  was  consistent  with  the  observationally-motivated  skeleton  model 

641  simulations  performed  by  Ogrosky  and  Stechmann  [2015].  Moreover,  the  leading  modes  of  the 

642  multivariate  EOF  analysis  from  Case-SST8N  had  greater  correlations  than  the  corresponding 

643  modes  from  Case-WP  while  also  containing  zonal  asymmetries  that  more  closely  resembled  the 

644  observed  Wheeler  and  Hendon  [2004]  RMM  EOFs.  Overall,  the  leading  EOFs  for  the 

645  multivariate  analysis  of  Case-SST8N  showed  a  wavenumber- 1  convectively-coupled  circulation 

646  anomaly  that  explained  15.7%  of  the  total  variance.  The  skeleton  multivariate  MJO  (SMM) 

647  index  also  relied  more  heavily  on  the  circulation  component,  as  is  the  case  in  observations  [e.g., 

648  Straub ,  2013],  given  the  lack  of  a  well-defined  univariate  EOF  mode  for  convective  heating  in 

649  the  model  data.  Correlation  coefficients  between  the  values  of  SMM1  and  SMM2  derived  from  a 

650  univariate  zonal  wind  analysis  (i.e.,  u  only)  and  the  full  SMM  index  (Ha  and  u)  were  in  excess 

65 1  of  0.98  when  considering  680  years  of  model  data  for  both  Case-WP  and  Case-SST8N. 

652  A  climatology  of  MJO  events  was  performed  using  multiple  forms  of  the  SMM  index  for 

653  the  skeleton  model.  It  was  shown  that  the  stochastic  skeleton  model  produces  a  number  of  MJO 

654  events  that,  for  such  a  simple  model,  compares  reasonably  well  with  the  observed  number  of 

655  MJO  events  documented  in  Stachnik  et  al.  [2015].  As  such,  the  average  MJO  event  length  in  the 

656  model  was  approximately  40  days  and  44  days  for  Case-WP  and  Case-SST8N,  respectively, 

657  which  is  only  slightly  higher  than  in  observations  (~35  days).  Distributions  of  the  average  event 

658  length  revealed  that  the  skeleton  model  does  a  good  job  of  simulating  the  shortest  40-50%  of 

659  MJO  events  and  that  the  average  length  is  skewed  by  the  occasional  presence  of  very  long  MJOs 

660  not  seen  in  nature,  suggesting  that  the  stochastic  parameters  used  in  the  model  are  well-tuned. 
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661  The  existence  of  additional  circumnavigating  and  long-lasting  MJOs  in  the  model  is  reasonable, 

662  however,  given  the  lack  of  damping  on  the  planetary  scale  and  the  omission  of  other  factors  such 

663  as  enhanced  friction  over  land  surfaces  that  may  help  tenninate  observed  events.  The  total 

664  number  of  MJOs  and  corresponding  event  lengths  were  relatively  insensitive  to  the  choice  of 

665  EOFs  (e.g.,  model  univariate,  multivariate,  or  the  observed  RMM  patterns)  used  to  calculate  the 

666  SMM  index. 

667  An  initial  attempt  was  made  to  quantity  the  geographic  distributions  of  MJO  event 

668  initiation  and  termination  based  on  the  large-scale  forcing.  Whereas  most  of  the  model  quantities 

669  were  relatively  stable  for  shorter  periods  that  mimicked  the  length  of  observations  (34  years),  the 

670  frequencies  of  MJO  initiation  in  a  given  RMM  phase  varied  strongly  over  shorter  time  periods 

67 1  and  required  a  much  longer  record  (680  years)  for  analysis.  While  the  simulated  climatology  for 

672  shorter  periods  may  occasionally  resemble  observations,  the  long-term  record  for  the  stochastic 

673  model  simulations  showed  nearly  equal  chance  for  an  MJO  event  to  initiate  or  terminate  in  any 

674  RMM  phase  despite  the  existence  of  zonal  variations  in  the  external  radiative  cooling  and 

675  moistening  terms  used  to  force  the  model.  The  deterministic  solutions  did  show  statistically 

676  significant  preferred  regions  of  MJO  initiation  and  termination  for  the  long-term  analysis,  though 

677  these  features  did  not  always  match  observations. 

678  The  previous  finding  presents  a  bit  of  conundrum,  considering  that  using  observed  SSTs 

679  to  force  the  model  results  in  a  more  realistic  representation  of  the  mean  convective  heating  and 

680  variance  that  is  in  agreement  with  Ogrosky  and  Stechmann  [2015]  while  also  producing 

68 1  multivariate  EOF  structures  that  more  closely  resemble  the  observed  RMM  patterns  taken  from 

682  Wheeler  and  Hendon  [2004].  Despite  this  fact,  the  stochastic  simulation  using  the  observed  SSTs 

683  provided  no  improvement  over  the  idealized  wann  pool  case  for  accurately  simulating  the 
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684  observed  geographic  variability  of  MJO  initiation  and  tennination  events.  The  effect  of 

685  stochasticity  here  appears  to  distribute  the  start  and  finish  of  individual  MJO  events  more  equally 

686  throughout  the  domain  while  maintaining  the  overall  amount  of  MJO  activity  in  the  model.  The 

687  MJO  events  initiated  and  tenninated  in  all  RMM  phases,  albeit  some  phases  were  preferred  for 

688  the  deterministic  model  simulations  that  did  not  necessarily  match  observations. 

689  While  the  basic  model  assumptions  and  observed,  annual  average  SSTs  may  be  sufficient 

690  to  generate  reasonable  looking  MJO  dynamics,  structures,  variability,  and  number  of  events, 

69 1  additional  factors  beyond  internal  tropical  variability  and  the  climatological  forcing  may  be 

692  required  to  produce  realistic  distributions  of  MJO  initiation  and  decay  events.  For  example, 

693  MJOs  may  initiate  from  any  number  of  mechanisms  including  influences  from  higher  latitudes 

694  such  as  cold  air  outbreaks  and  Rossby  wave  trains  [e.g.,  Hsu  et  al.,  1990;  Matthews  and  Kiladis, 

695  1999;  Ray  and  Zhang,  2010;  Wang  et  al.,  2012]  or  topographic  effects  [e.g.,  Hsu  and  Lee,  2005], 

696  all  of  which  have  geographic  variability  that  are  not  accounted  for  in  the  current  version  of  the 

697  skeleton  model.  Future  work  will  continue  to  evaluate  long-term  records  of  simulated  MJO 

698  behavior  from  the  skeleton  model  as  a  function  of  the  shape  and  strength  of  the  tropical  warm 

699  pool  in  order  to  determine  if  the  model  is  capable  of  identifying  shifts  in  the  location  or  lifetimes 

700  of  MJO  events  as  a  function  of  ENSO  [e.g.,  Pohl  and  Matthews,  2007]  or  Indian  Ocean  dipole 

701  (IOD)SSTs. 

702  Including  stochastic  effect  within  the  convective  parameterization  were  also  necessary  for 

703  the  model  to  produce  more  realistic  estimates  of  the  time  scales  for  MJO  event  initiation  and 

704  tennination.  The  introduction  of  the  stochastic  parameterization  better  accounts  for  the  random 

705  component  of  the  growth  and  decay  of  unresolved  synoptic  and  mesoscale  convective  processes 

706  that  occur  within  the  MJO  envelope.  These  processes  are  needed  here,  and  likely  in  nature,  to 
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707  accelerate  the  time  scales  of  initiation  and  decay  that  are  otherwise  too  slow  through  the 

708  neutrally  stable  wave  interactions  and  large-scale  moisture  depletion  that  occurs  in  the 

709  detenninistic  version  of  the  skeleton  model.  In  fact,  almost  every  one  of  the  stochastic  model 

710  simulations  produced  a  more  realistic  looking  MJO  initiation  and  tennination  composite  with 

711  smaller  RMSE  values  regardless  of  the  choice  of  EOFs  when  calculating  the  principal 

712  components  or  background  forcing  when  compared  to  its  corresponding  detenninistic  solution 

713  (not  shown). 

714  Finally,  we  examined  the  large-scale  conditions  associated  with  the  termination  of  MJO 

715  events  in  the  skeleton  model  over  the  RMM  phase  6  domain.  As  previously  reported  through  a 

716  visual  analysis  of  MJO  time  series  in  Majda  and  Stechmann  [2009,  2011]  and  Thual  at  al. 

Ill  [2014],  the  skeleton  model  produced  MJO  events  that  either  stalled  or  experienced  decay 

718  immediately  upstream  or  directly  over  the  warm  pool  maximum  as  observed  in  nature.  However, 

719  it  is  important  to  note  that  these  successes  may  not  be  for  the  correct  physical  reasons,  as  the 

720  observed  MJO  encounters  a  number  of  additional  factors  that  may  weaken  it  over  this  domain 

72 1  that  are  not  accounted  for  in  the  skeleton  model  such  as  stronger  diurnal  cycles,  topographic 

722  effects,  and  enhanced  boundary  layer  friction  over  land.  Moreover,  our  analysis  of  the  precursor 

723  conditions  associated  with  MJO  decay  in  the  skeleton  model  indicated  that  higher  wavenumber 

724  ( k  ~  12-13)  oscillations  embedded  within  the  large-scale  flow  could  be  important  for  MJO 

725  tennination.  Whether  the  development  of  these  waves  are  either  suppressed  by  continuing  MJO 

726  events,  or  alternatively,  the  emergence  of  eastward  propagating  waves  contribute  towards  MJO 

727  event  termination  remains  unknown.  Our  future  work  plans  to  examine  the  wave  power  spectra 

728  associated  with  MJO  tennination  events  in  observations  and  the  skeleton  model  to  better 

729  elucidate  what  mechanisms  may  be  responsible  for  MJO  initiation  and  decay. 
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889  Figure  Captions 

890  Figure  1.  The  shape  of  (a)  the  idealized  warm  pool  (i.e.,  Case-WP)  and  (b)  annual  average 

891  (1971-2000)  NOAA  01  SSTs  v2  used  to  force  the  skeleton  model.  The  observed  SSTs  have  been 

892  meridionally  averaged  from  10°S-10°N  with  a  Poisson  relaxation  and  interpolated  over  missing 

893  data  (i.e.,  land).  The  solid  line  in  (b)  indicates  the  average  SSTs  from  all  waves  downscaled  to 

894  the  model  resolution  of  625  km.  The  dashed  profile  in  (b)  has  been  smoothed  to  only  retain  the 

895  first  eight  Fourier  modes  (Case-SST8N). 

896 

897  Figure  2.  (a)-(c)  Hovmoller  diagrams  of  the  anomalous  20-100  day  filtered  skeleton  model 

898  variables  for  a  600  day  period  taken  from  the  idealized  warm  pool  simulation  (Case-WP).  All 

899  variables  have  daily  resolution  and  are  shown  in  dimensional  units  for  the  reconstructed  (a) 

900  convective  heating  (Ha),  (b)  zonal  wind  (w),  and  (c)  moisture  (q)  variables,  (d)-(f)  As  in  (a)-(c), 

901  but  for  the  model  simulation  using  the  truncated,  observed  SSTs  (Case-SST8N). 

902 

903  Figure  3.  Zonal  wavenumber-frequency  power  spectra  of  the  (a)  convective  heating  (Ha,  K  day' 

904  ')  and  (b)  zonal  wind  (u,  m  s"1)  for  the  idealized  warm  pool  case  (Case-WP)  from  the  skeleton 

905  model.  Shading  indicates  the  base- 10  logarithm  for  the  reconstructed  and  dimensional  variables 

906  taken  at  the  equator  for  approximately  34  years  of  model  data.  The  black  circles  mark  the 

907  dispersion  curves  from  a  separate  stability  analysis  of  the  linear  waves  in  the  skeleton  model  in 

908  the  Indian  Ocean  (67.5°E).  The  dashed,  black  lines  indicate  periods  of  30  and  90  days,  (c)-(d)  As 

909  in  (a)-(b),  but  for  the  model  simulation  using  the  truncated,  observed  SSTs  (Case-SST8N). 

910 
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911  Figure  4.  (a)-(b)  The  first  two  leading  modes  from  a  univariate  EOF  analysis  of  the  20-100  day 

912  filtered  convective  heating  (Ha)  anomalies  along  the  equator  from  approximately  340  years  of 

913  the  idealized  warm  pool  simulation  (Case-WP).  (c)-(d)  As  in  (a)-(b),  but  for  a  univariate  EOF 

914  analysis  of  the  zonal  wind  (u).  The  percent  variance  explained  by  each  mode  is  shown  in  the 

9 1 5  upper-right  of  each  panel. 

916 

917  Figure  5.  (a)-(b)  The  first  two  leading  modes  from  a  multivariate  EOF  analysis  of  the  20-100 

918  day  filtered  convective  heating  (—Ha)  and  zonal  wind  (u)  model  anomalies  along  the  equator 

919  from  approximately  340  years  of  the  idealized  warm  pool  simulation  (Case-WP).  The  order  of 

920  EOF1  and  EOF2  have  been  interchanged  in  addition  to  switching  the  sign  of  the  convective 

921  heating  term  to  better  match  the  Wheeler  and  Hendon  [2004]  EOFs.  (c)-(d)  The  first  leading  two 

922  modes  from  the  observed  multivariate  analysis  of  Wheeler  and  Hendon  [2004]  showing  the 

923  15°S-15°N  EOF  structures  of  OLR  and  zonal  wind  at  850  hPa  and  200  hPa  (u850  and  u200, 

924  respectively).  The  percent  variance  explained  by  each  mode  is  shown  in  the  upper-right  of  each 

925  panel. 

926 

927  Figure  6.  (a)-(b)  As  in  Figure  5,  but  for  the  model  simulation  using  the  truncated,  observed  SSTs 

928  (Case-SST8N).  Panels  (c)-(d)  are  identical  to  Figure  5  to  better  facilitate  comparison  between  the 

929  modeled  and  observed  EOFs. 

930 

93 1  Figure  7.  Fag  correlations  between  the  corresponding  PCI  and  PC2  from  the  multivariate  EOF 

932  analysis  of  approximately  340  years  of  the  model  simulation  using  the  truncated,  observed  SSTs 

933  (Case-SST8N). 
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934  Figure  8.  Examples  of  two  MJO  events  using  the  skeleton  multivariate  MJO  (SMM)  index  and 

935  the  event  criteria  defined  in  Stachnik  el  al.  [2015]  for  the  model  simulation  using  an  idealized 

936  warm  pool  (Case-WP).  The  daily  values  of  SMM1  and  SMM2  (i.e.,  the  principal  components) 

937  are  plotted  as  small  dots  for  (a)  day  30561-30619  and  (b)  day  32662-32734  in  the  2-dimensional 

938  phase  space  obtained  when  projecting  the  model  data  onto  the  observed  EOF  patterns  of  Wheeler 

939  and  Hendon  [2004].  The  initiation  of  a  primary  MJO  event  is  shown  in  gray  and  indicated  with  a 

940  large,  circumscribed  X  while  the  corresponding  termination  event  is  shown  in  red  for  each  panel. 

941  Continuing  MJO  events  are  shown  in  blue,  green,  and  orange  in  (a),  with  the  subsequent  history 

942  of  the  SMM  index  highlighted  through  the  following  four  phases.  The  beginning  of  each 

943  continuing  event  are  shown  with  large,  blue  markers  in  (b)  along  with  the  start  of  a 

944  circumnavigating  event  and  associated  path  highlighted  in  gold. 

945 

946  Figure  9.  Cumulative  probability  distribution  functions  of  the  observed  and  modeled  MJO  event 

947  lengths  from  the  stochastic  model  simulation  using  the  truncated,  observed  SSTs  (Case-SST8N), 

948  and  a  deterministic  model  simulation  using  the  same  SSTs  (Case-SST8Ndet).  The  observed 

949  distribution  is  based  on  the  historical  Wheeler  and  Hendon  [2004]  RMM  index  from  1979-2012 

950  and  is  represented  by  a  thick,  black  line.  The  corresponding  model  event  length  distributions  for 

95 1  similar  34  year  periods  are  shown  using  thin,  blue  and  red  lines  for  the  stochastic  and 

952  deterministic  solutions,  respectively. 

953 

954  Figure  10.  Radial-frequency  plots  showing  the  percent  occurrence  of  MJO  events  as  a  function 

955  of  RMM  phase  for  (a)  observations  (1979-2012)  and  single  34  year  periods  from  the  stochastic 

956  model  simulations  using  (b)  an  idealized  warm  pool  (Case-WP)  and  (c)  truncated,  observed  SSTs 
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957  (Case-SST8N).  (d)  As  in  (c),  but  for  the  deterministic  model  simulation  using  the  truncated, 

958  observed  SSTs  (Case-SST8Ndet).  The  black  dots  plotted  at  the  center  angle  of  each  phase 

959  indicate  the  frequency  of  occurrence  and  are  connected  only  for  the  ease  of  viewing  primary, 

960  continuing,  and  terminal  MJO  events  (gray,  blue,  and  red  lines,  respectively).  Range  rings 

961  represent  the  1%,  5%,  10%,  15%,  and  20%  frequencies  of  occurrence  from  the  center  to  edge  of 

962  the  diagram,  respectively. 

963 

964  Figure  11.  The  mean  frequency  and  uncertainty  estimates  of  primary  (gray),  continuing  (blue), 

965  and  terminal  (red)  MJO  events  for  each  RMM  phase  in  (a)  observations  and  (b)-(f) 

966  approximately  680  years  of  model  data.  Simulations  using  an  idealized  warm  pool  and  the 

967  truncated,  observed  SSTs  (Case-WP  and  Case-SST8N,  respectively)  are  shown  in  panels  (b)-(c) 

968  with  a  special  simulation  using  a  homogeneous  background  state  (Case-HG)  in  (d).  Results  from 

969  a  deterministic  version  of  the  skeleton  model  using  the  idealized  warm  pool  (Case-WPdet)  and 

970  truncated,  observed  SSTs  (Case-SST8Ndet)  are  shown  in  (e)-(f).  The  dashed  line  indicates  the 

971  12.5%  frequency  and  represents  the  equal  likelihood  of  an  MJO  event  occurring  in  any  phase. 

972  Error  bars  indicate  the  95%  uncertainty  estimates  and  are  calculated  from  a  binomial  proportion 

973  confidence  interval  (see  text  for  details).  The  total  number  of  primary/terminal  and  continuing 

974  events  are  indicated  in  the  top  right  of  each  panel. 

975 

976  Figure  12.  Initiation  composites  of  the  RMM1-RMM2  paths  for  (a)  observations  (1979-2012) 

977  and  SMM1-SMM2  for  an  equivalent  34  year  period  from  the  stochastic  model  simulations  using 

978  (b)  an  idealized  warm  pool  (Case-WP)  and  (c)  the  truncated,  observed  SSTs  (Case-SST8N).  (d) 

979  As  in  (c),  but  for  the  deterministic  model  simulation  using  the  truncated  and  observed  SSTs 


45 


980  (Case-SST8Ndet).  The  composite  values  for  a  primary  event  are  marked  with  a  circumscribed  X 

981  for  each  phase.  The  composite  history  before  (3  days)  and  after  (12  days)  the  initiation  event 

982  (i.e.,  first  day  after  the  combined  amplitude  grows  above  unity)  are  shown  in  different  colors  for 

983  each  phase.  The  composite  daily  RMM  and  SMM  values  are  plotted  using  small  dots,  while 

984  closed  squares  appear  every  three  days  for  reference. 

985 

986  Figure  13.  As  in  Figure  12,  but  for  MJO  termination  events.  The  composite  values  for  a 

987  termination  event  are  marked  with  a  circumscribed  X  for  each  phase.  The  composite  history 

988  before  (12  days)  and  after  (3  days)  the  termination  event  (i.e.,  last  day  before  the  combined 

989  amplitude  decays  below  unity)  are  shown  in  different  colors  for  each  phase. 

990 

991  Figure  14.  Longitude-lag  composites  of  the  20-100  day  filtered  anomalous  (a)-(b)  convective 

992  heating  (Ha),  (d)-(e)  moisture  variable  (q),  and  (g)-(h)  zonal  wind  (u)  from  the  model  simulation 

993  using  the  truncated,  observed  SSTs  (Case-SST8N)  for  RMM  phase  6  continuing  and  terminal 

994  MJO  events.  All  data  are  taken  along  the  equator  and  composited  for  approximately  340  years. 

995  Differences  between  terminal  and  continuing  MJO  events  are  shown  in  panels  (c),  (f),  and  (i)  for 

996  Ha,  q,  and  u,  respectively.  The  stippling  in  the  difference  panels  indicates  the  presence  of  a 

997  statistically  significant  different  local  mean  for  the  continuing  and  terminating  events  (see 

998  Stachnik  el  al.  [2015]  for  further  details).  The  number  of  events  used  for  each  composite  (N)  is 

999  listed  in  the  top-right  corner  of  each  panel. 

1000 
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1001  Figure  15.  As  in  Figure  14,  but  for  the  (a)-(c)  zonal  convergence  and  (d)-(f)  potential 

1002  temperature  ( 0).  (g)-(i)  As  in  (d)-(f),  but  for  a  shorter,  single  34  year  period  of  skeleton  model 

1003  data. 

1004 

1005  Figure  16.  Box- whisker  diagrams  showing  the  distribution  of  (a)  zonal  wind  (u)  and  (b) 

1006  moisture  variable  ( q )  anomalies  for  RMM  phase  6  continuing  and  terminating  MJO  events 

1007  averaged  over  the  local  domain  (120-160°E)  at  day  -10.  Box  lines  indicate  the  25%,  50%,  and 

1008  75%  percentiles  with  whiskers  indicating  the  minimum  and  maximum  values  for  all  events.  The 

1009  blue  circles  and  red  triangles  represent  the  local-domain  average  values  for  individual  continuing 

1010  and  terminating  MJO  events,  respectively.  The  mean  quantities  for  each  event  type  are  overlaid 

1011  with  a  black  star.  All  data  are  taken  along  the  equator  and  represent  34  years  of  simulated  MJO 

1012  events  using  the  truncated,  observed  SSTs  (Case-SST8N). 

1013 
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1014  Table  1.  Maximum  correlation  coefficient  for  the  longitudinal  structure  of  the  leading  two  EOFs 

1015  from  all  lags  between  the  skeleton  model  and  observed  Wheeler  and  Hendon  [2004]  RMM  index 


EOF  Variable 

EOF 

Univariate  (u) 

Multivariate  (w) 

Multivariate  (—Ha) 

Case-WP 

EOF1 

0.85 

0.83 

0.70 

EOF2 

0.74 

0.72 

0.68 

Case-SST8N 

EOF1 

0.98 

0.98 

0.89 

EOF2 

0.93 

0.94 

0.80 

1016 
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1017  Table  2.  Number  of  MJO  events  and  average  event  length  from  observations  (1979-2012) a  and 

1018  an  equivalent  34  year  periodb  from  the  stochastic  skeleton  model  for  the  simulation  using  an 

1019  idealized  warm  pool  (Case-WP)  as  a  function  of  EOFs  used  to  calculate  the  SMM  index 


Event  type  (#) 

EOFs 

Primary 

Continuing 

Circum¬ 

navigating 

Terminal 

Event 
length  (d) 

Univariate  (u) 
Multivariate 

89.3±4.3 

309.2±23.4 

21.2±3.1 

89.3±4.3 

39.5+1.6 

(Model) 

Multivariate 

91.4±3.9 

341±22.9 

23.5±3.0 

91.4±3.9 

39.9±1.4 

(WH04C) 

97.7±3.5 

343.4±17.9 

23.7±2.0 

97.7±3.5 

40.2±1.1 

Observations 

154 

330 

15 

154 

34.9 

1020 

1021  aThe  observed  MJO  event  numbers  are  taken  from  Stachnik  el  al.  [2015] 

1022  bThe  model  statistics  and  uncertainty  estimates  are  generated  from  680  years  of  model  data 

1023  comprising  20  independent  samples  of  34  year  periods 

1024  CWH04  refers  to  the  case  in  which  the  observed  Wheeler  and  Hendon  [2004]  EOFs  are  used  in 

1025  place  of  the  model  quantities  to  calculate  the  SMM  index 
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1026  Table  3.  As  in  Table  2,  but  for  the  simulation  using  the  truncated,  observed  SSTs  (Case-SST8N) 


Event  type  (#) 

EOFs 

Primary 

Continuing 

Circum¬ 

navigating 

Terminal 

Event 
length  (d) 

Univariate  (u) 
Multivariate 

98.3±2.8 

443.5±31.5 

33.4±3.6 

98.3±2.8 

43.3±1.7 

(Model) 

Multivariate 

98.4±3.2 

466.8±31.2 

35.5±3.4 

98.4±3.2 

43.8±1.5 

(WH04) 

97.7±3.6 

477.9±33.2 

36.6±3.5 

97.7±3.6 

44.2±1.6 

Observations 

154 

330 

15 

154 

34.9 

1027 
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1028  Table  4.  Number  of  MJO  events  and  average  event  length  from  observations  (1979-2012)  and  an 


1029  equivalent  34  year  period  from  the  stochastic  and  deterministic  skeleton  model  simulations  using 

1030  the  idealized  warm  pool  (Case-WP  and  Case-WPdet)  and  truncated,  observed  SSTs  (Case- 
103  1  SST8N  and  Case-SST8Ndet)  projected  on  the  observed  Wheeler  and  Hendon  [2004]  RMM 
1032  EOFs 


Dataset 

Event  type  (#) 

Primary 

Continuing 

Circum¬ 

navigating 

Terminal 

Event 
length  (d) 

Case-WP 

97.7±3.5 

343.4±17.9 

23.7±2.0 

97.7+3.5 

40.2+1.1 

Case-SST8N 

97.7±3.6 

477.9±33.2 

36.6±3.5 

97.7+3.6 

44.2+1.6 

Case-WPdet 

90.8±2.3 

797.3±29.0 

75.3±3.7 

90.8+2.3 

61.5+2.2 

Case-SST8Ndet 

87.4±1.9 

1316.7±26.6 

140.1+3.2 

87.4+1.9 

77.7+1.8 

Observations 

154 

330 

15 

154 

34.9 

1033 
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1038 
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Figure  1.  The  shape  of  (a)  the  idealized  warm  pool  (i.e.,  Case-WP)  and  (b)  annual  average 
(1971-2000)  NOAA  OI  SSTs  v2  used  to  force  the  skeleton  model.  The  observed  SSTs  have  been 
meridionally  averaged  from  10°S-10°N  with  a  Poisson  relaxation  and  interpolated  over  missing 
data  (i.e.,  land).  The  solid  line  in  (b)  indicates  the  average  SSTs  from  all  waves  downscaled  to 
the  model  resolution  of  625  km.  The  dashed  profile  in  (b)  has  been  smoothed  to  only  retain  the 
first  eight  Fourier  modes  (Case-SST8N). 
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Figure  2.  (a)-(c)  Hovmoller  diagrams  of  the  anomalous  20-100  day  filtered  skeleton  model 
variables  for  a  600  day  period  taken  from  the  idealized  warm  pool  simulation  (Case-WP).  All 
variables  have  daily  resolution  and  are  shown  in  dimensional  units  for  the  reconstructed  (a) 
convective  heating  (Ha),  (b)  zonal  wind  ( u ),  and  (c)  moisture  ( q )  variables,  (d)-(f)  As  in  (a)-(c), 
but  for  the  model  simulation  using  the  truncated,  observed  SSTs  (Case-SST8N). 
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Figure  3.  Zonal  wavenumber-frequency  power  spectra  of  the  (a)  convective  heating  (Ha,  K  day" 
')  and  (b)  zonal  wind  (u,  m  s'1)  for  the  idealized  warm  pool  case  (Case-WP)  from  the  skeleton 
model.  Shading  indicates  the  base- 10  logarithm  for  the  reconstructed  and  dimensional  variables 
taken  at  the  equator  for  approximately  34  years  of  model  data.  The  black  circles  mark  the 
dispersion  curves  from  a  separate  stability  analysis  of  the  linear  waves  in  the  skeleton  model  in 
the  Indian  Ocean  (67.5°E).  The  dashed,  black  lines  indicate  periods  of  30  and  90  days,  (c)-(d)  As 
in  (a)-(b),  but  for  the  model  simulation  using  the  truncated,  observed  SSTs  (Case-SST8N). 
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1060  Figure  4.  (a)-(b)  The  first  two  leading  modes  from  a  univariate  EOF  analysis  of  the  20-100  day 

1061  filtered  convective  heating  (Ha)  anomalies  along  the  equator  from  approximately  340  years  of 

1062  the  idealized  warm  pool  simulation  (Case-WP).  (c)-(d)  As  in  (a)-(b),  but  for  a  univariate  EOF 

1063  analysis  of  the  zonal  wind  ( u ).  The  percent  variance  explained  by  each  mode  is  shown  in  the 

1064  upper-right  of  each  panel. 
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Figure  5.  (a)-(b)  The  first  two  leading  modes  from  a  multivariate  EOF  analysis  of  the  20-100 
day  filtered  convective  heating  (—Ha)  and  zonal  wind  ( u )  model  anomalies  along  the  equator 
from  approximately  340  years  of  the  idealized  warm  pool  simulation  (Case-WP).  The  order  of 
EOF1  and  EOF2  have  been  interchanged  in  addition  to  switching  the  sign  of  the  convective 
heating  term  to  better  match  the  Wheeler  and  Hendon  [2004]  EOFs.  (c)-(d)  The  first  leading  two 
modes  from  the  observed  multivariate  analysis  of  Wheeler  and  Hendon  [2004]  showing  the 
15°S-15°N  EOF  structures  of  OLR  and  zonal  wind  at  850  hPa  and  200  hPa  (u850  and  u200, 
respectively).  The  percent  variance  explained  by  each  mode  is  shown  in  the  upper-right  of  each 
panel. 
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Figure  6.  (a)-(b)  As  in  Figure  5,  but  for  the  model  simulation  using  the  truncated,  observed  SSTs 
(Case-SST8N).  Panels  (c)-(d)  are  identical  to  Figure  5  to  better  facilitate  comparison  between  the 
modeled  and  observed  EOFs. 
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1083  Figure  7.  Lag  correlations  between  the  corresponding  PCI  and  PC2  from  the  multivariate  EOF 

1084  analysis  of  approximately  340  years  of  the  model  simulation  using  the  truncated,  observed  SSTs 

1085  (Case-SST8N). 
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(a)  SMM  Index  for  Day  30561-3061 9 
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Phase  7  (Western  Pacific)  Phase  6 
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Figure  8.  Examples  of  two  MJO  events  using  the  skeleton  multivariate  MJO  (SMM)  index  and 
the  event  criteria  defined  in  Stachnik  et  al.  [2015]  for  the  model  simulation  using  an  idealized 
warm  pool  (Case-WP).  The  daily  values  of  SMM1  and  SMM2  (i.e.,  the  principal  components) 
are  plotted  as  small  dots  for  (a)  day  30561-30619  and  (b)  day  32662-32734  in  the  2-dimensional 
phase  space  obtained  when  projecting  the  model  data  onto  the  observed  EOF  patterns  of  Wheeler 
and  Hendon  [2004].  The  initiation  of  a  primary  MJO  event  is  shown  in  gray  and  indicated  with  a 
large,  circumscribed  X  while  the  corresponding  termination  event  is  shown  in  red  for  each  panel. 
Continuing  MJO  events  are  shown  in  blue,  green,  and  orange  in  (a),  with  the  subsequent  history 
of  the  SMM  index  highlighted  through  the  following  four  phases.  The  beginning  of  each 
continuing  event  are  shown  with  large,  blue  markers  in  (b)  along  with  the  start  of  a 
circumnavigating  event  and  associated  path  highlighted  in  gold. 
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1100 

1101  Figure  9.  Cumulative  probability  distribution  functions  of  the  observed  and  modeled  MJO  event 

1 102  lengths  from  the  stochastic  model  simulation  using  the  truncated,  observed  SSTs  (Case-SST8N), 

1103  and  a  deterministic  model  simulation  using  the  same  SSTs  (Case-SST8Ndet).  The  observed 

1 104  distribution  is  based  on  the  historical  Wheeler  and  Hendon  [2004]  RMM  index  from  1979-2012 

1105  and  is  represented  by  a  thick,  black  line.  The  corresponding  model  event  length  distributions  for 

1106  similar  34  year  periods  are  shown  using  thin,  blue  and  red  lines  for  the  stochastic  and 

1107  deterministic  solutions,  respectively. 
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Figure  10.  Radial-frequency  plots  showing  the  percent  occurrence  of  MJO  events  as  a  function 
of  RMM  phase  for  (a)  observations  (1979-2012)  and  single  34  year  periods  from  the  stochastic 
model  simulations  using  (b)  an  idealized  warm  pool  (Case-WP)  and  (c)  truncated,  observed  SSTs 
(Case-SST8N).  (d)  As  in  (c),  but  for  the  deterministic  model  simulation  using  the  truncated, 
observed  SSTs  (Case-SST8Ndet).  The  black  dots  plotted  at  the  center  angle  of  each  phase 
indicate  the  frequency  of  occurrence  and  are  connected  only  for  the  ease  of  viewing  primary, 
continuing,  and  terminal  MJO  events  (gray,  blue,  and  red  lines,  respectively).  Range  rings 
represent  the  1%,  5%,  10%,  15%,  and  20%  frequencies  of  occurrence  from  the  center  to  edge  of 
the  diagram,  respectively. 
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Figure  11.  The  mean  frequency  and  uncertainty  estimates  of  primary  (gray),  continuing  (blue), 
and  terminal  (red)  MJO  events  for  each  RMM  phase  in  (a)  observations  and  (b)-(f) 
approximately  680  years  of  model  data.  Simulations  using  an  idealized  warm  pool  and  the 
truncated,  observed  SSTs  (Case-WP  and  Case-SST8N,  respectively)  are  shown  in  panels  (b)-(c) 
with  a  special  simulation  using  a  homogeneous  background  state  (Case-HG)  in  (d).  Results  from 
a  deterministic  version  of  the  skeleton  model  using  the  idealized  warm  pool  (Case-WPdet)  and 
truncated,  observed  SSTs  (Case-SST8Ndet)  are  shown  in  (e)-(f).  The  dashed  line  indicates  the 
12.5%  frequency  and  represents  the  equal  likelihood  of  an  MJO  event  occurring  in  any  phase. 
Error  bars  indicate  the  95%  uncertainty  estimates  and  are  calculated  from  a  binomial  proportion 
confidence  interval  (see  text  for  details).  The  total  number  of  primary/terminal  and  continuing 
events  are  indicated  in  the  top  right  of  each  panel. 
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Figure  12.  Initiation  composites  of  the  RMM1-RMM2  paths  for  (a)  observations  (1979-2012) 
and  SMM1-SMM2  for  an  equivalent  34  year  period  from  the  stochastic  model  simulations  using 
(b)  an  idealized  warm  pool  (Case-WP)  and  (c)  the  truncated,  observed  SSTs  (Case-SST8N).  (d) 
As  in  (c),  but  for  the  deterministic  model  simulation  using  the  truncated  and  observed  SSTs 
(Case-SST8Ndet).  The  composite  values  for  a  primary  event  are  marked  with  a  circumscribed  X 
for  each  phase.  The  composite  history  before  (3  days)  and  after  (12  days)  the  initiation  event 
(i.e.,  first  day  after  the  combined  amplitude  grows  above  unity)  are  shown  in  different  colors  for 
each  phase.  The  composite  daily  RMM  and  SMM  values  are  plotted  using  small  dots,  while 
closed  squares  appear  every  three  days  for  reference. 
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Figure  13.  As  in  Figure  12,  but  for  MJO  termination  events.  The  composite  values  for  a 
termination  event  are  marked  with  a  circumscribed  X  for  each  phase.  The  composite  history 
before  (12  days)  and  after  (3  days)  the  termination  event  (i.e.,  last  day  before  the  combined 
amplitude  decays  below  unity)  are  shown  in  different  colors  for  each  phase. 
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Figure  14.  Longitude-lag  composites  of  the  20-100  day  filtered  anomalous  (a)-(b)  convective 
heating  (Ha),  (d)-(e)  moisture  variable  ( q ),  and  (g)-(h)  zonal  wind  ( u )  from  the  model  simulation 
using  the  truncated,  observed  SSTs  (Case-SST8N)  for  RMM  phase  6  continuing  and  terminal 
MJO  events.  All  data  are  taken  along  the  equator  and  composited  for  approximately  340  years. 
Differences  between  terminal  and  continuing  MJO  events  are  shown  in  panels  (c),  (f),  and  (i)  for 
Ha,  q,  and  u,  respectively.  The  stippling  in  the  difference  panels  indicates  the  presence  of  a 
statistically  significant  different  local  mean  for  the  continuing  and  terminating  events  (see 
Stachnik  et  al.  [2015]  for  further  details).  The  number  of  events  used  for  each  composite  (N)  is 
listed  in  the  top-right  corner  of  each  panel. 
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Figure  15.  As  in  Figure  14,  but  for  the  (a)-(c)  zonal  convergence  and  (d)-(f)  potential 
temperature  ( 0).  (g)-(i)  As  in  (d)-(f),  but  for  a  shorter,  single  34  year  period  of  skeleton  model 
data. 
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Figure  16.  Box-whisker  diagrams  showing  the  distribution  of  (a)  zonal  wind  (u)  and  (b) 
moisture  variable  ( q )  anomalies  for  RMM  phase  6  continuing  and  terminating  MJO  events 
averaged  over  the  local  domain  (120-160°E)  at  day  -10.  Box  lines  indicate  the  25%,  50%,  and 
75%  percentiles  with  whiskers  indicating  the  minimum  and  maximum  values  for  all  events.  The 
blue  circles  and  red  triangles  represent  the  local-domain  average  values  for  individual  continuing 
and  terminating  MJO  events,  respectively.  The  mean  quantities  for  each  event  type  are  overlaid 
with  a  black  star.  All  data  are  taken  along  the  equator  and  represent  34  years  of  simulated  MJO 
events  using  the  truncated,  observed  SSTs  (Case-SST8N). 
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